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Abstract 



The environment in which bremsstrahlung and pair creation occurs can strongly 
affect cross sections for these processes. Because ultra-relativistic electromagnetic 
interactions involve very small longitudinal momentum transfers, the reactions occur 
gradually, spread over long distances. During this time, even relatively weak factors 
can accumulate enough to disrupt the interaction. 

This review will discuss a variety of factors which can suppress bremsstrahlung 
and pair production, as well as related effects involving beamstrahlung and QCD 
processes. After surveying different theoretical approaches, experimental measure- 
ments will be covered. Recent accurate measurements by the SLAC E-146 collab- 
oration will be highlighted, along with several recent theoretical works relating to 
the experiment. 

(To appear in Reviews of Modern Physics) 
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I. INTRODUCTION 

Bremsstrahlung and pair creation are two of the most common high energy electromag- 
netic processes, and the interaction cross sections are well known (Bethe and Heitler, 1934). 
However, it is much less well known that these cross sections can change dramatically de- 
pending on the environment in which the interaction occurs. The SLAC E-146 collaboration 
observed photon intensities as low as 1/4 of that predicted by Bethe and Heitler, and much 
larger suppression is possible. 
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The cross sections change because the kinematics of the processes dictates the interac- 
tions are spread over a significant distance, in contrast to the conventional picture where 
reactions occur at a single point. The momentum transfer between a highly relativistic 
interacting particle and the target nucleus can be small, especially along the direction of 
particle motion. When this longitudinal momentum transfer is small, the uncertainty prin- 
ciple dictates that the interaction is spread out over a distance, known as the formation 
length for particle production, or, more generally, as the coherence length. 

If the medium in the neighborhood of the interaction has enough influence on the inter- 
acting particle during its passage through the formation zone, then the pair production or 
bremsstrahlung can be suppressed. When the formation length is long, weak but cumulative 
factors can be important. Some of the factors that may suppress electromagnetic radiation 
are multiple scattering, photon interactions with the medium (coherent forward Compton 
scattering), and magnetic fields. In crystals, where the atoms are arranged in ordered rows, 
a large variety of effects can suppress or enhance radiation. 

The formation length has a number of interesting physical interpretations. It is the 
wavelength of the exchanged virtual photon. It is also the distance required for the final 
state particles to separate enough (by a Compton wavelength) that they act as separate 
particles. It is also the distance over which the amplitudes from several interactions can add 
coherently to the total cross section. All of these pictures can be helpful in understanding 
different aspects of suppression mechanisms. 

The formation length was first considered by Ter-Mikaelian in 1952 (Feinberg, 1994). It 
occurs in both classical and quantum mechanical calculations. Classically, all electromag- 
netic processes have a phase factor 

$ = exp {i[iut - k ■ r(t)]} (1) 

where ui is the photon frequency, t is the time, k is the photon wave vector, and r(t) is the 
electron position. The formation length is the distance over which this phase factor maintains 
coherence; the electron and photon have slightly different velocities, which causes them to 
gradually separate. Coherence is maintained over a distance for which ut — k ■ r(t) ~ 1. 
Around an isolated atom u = |k|c, so, neglecting the (small) angle between k and r, for a 
straight trajectory 

where E is the electron energy, m is the electron mass, and c the speed of light. The fO 
subscript denotes the unsuppressed (free space) formation length. Electromagnetic inter- 
actions are spread over this distance. If E is large and u small, lf can be very long. For 
example, for a 25 GeV electron radiating a 10 MeV bremsstrahlung photon, lf = 100/im. 
This distance was the basis for the suppression calculations by Landau and Pomeranchuk, 
discussed in Sec. II. 

The formation length also occurs in quantum mechanical calculations. Electron 
bremsstrahlung, where an electron interacts with a nucleus and emits a real photon is a 
good example. As the electron accelerates, part of its surrounding virtual photon field 
shakes loose. Neglecting the photon emission angle and electron scattering, the momentum 
transfer in the longitudinal (+z) direction is 
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?|| =Pe-p' e - P ,= y/(E/c)* -"- (mc) 2 - yj ((E — k)/c) 2 - (mc) 2 - fc/c (3) 



where p e and p' e are the electron momenta before and after the interaction respectively 
and the photon momentum p 7 = k/c, with k = fiuj = h\k\/c the photon energy. For 
7 = E/mc 2 3> 1 this simplifies to 

m 2 c 3 k . . 

Here, we have used the small y approximation y/\ — y ~ l — y/2. For high energy electrons 
emitting low energy photons, q» can become very small. For the above example, q» = 
0.002 eV/c. Because q\\ is so small, the uncertainty principle requires that the radiation 
must take place over a long distance: 

l h _ 2hE(E - k) 
^° 5|| m 2 c 3 k 

For k <^ E, this matches Eq. (2); the momentum transfer approach is needed only for 
k ~ E. 

This formation length appears in most electromagnetic processes including pair pro- 
duction, transition radiation, Cerenkov radiation and synchrotron radiation. Usually, q\\ 
acquires additional terms to account for specific interactions with the target medium. This 
review will discuss a number of these additions, and show how they affect bremsstrahlung 
and pair production. Formation lengths also apply to processes involving other forces; this 
review will consider a few non-electromagnetic interactions, showing how kinematics can 
dictate that diverse reactions share many fundamental traits. 

The formation length is the distance over which the interaction amplitudes can add 
coherently. So, if an external force doubles q\\, thereby halving If, then the radiation intensity 
is also halved, from ~ I 2 to ~ 2 • (lf/2) 2 as the trajectory splits into two independent 
emitting regions. If an electron traversing a distance D radiates, the trajectory acts as 
D/lf independent emitters, each with a strength proportional to \lf\ 2 , so the total radiation 
is proportional to If. As external factors increase q\\ and reduce //, the radiation drops 
proportionally. 

A similar coherence distance h/q± limits the perpendicular distance over which coherent 
addition is possible. However, because q± ^> q\\, this distance is much smaller, and of lesser 
interest here. 

This article begins by considering separate classical and quantum mechanical (semi- 
classical) calculations of suppression due to multiple scattering. Suppression due to photon 
interactions with the medium, pair creation and external magnetic fields will then be simi- 
larly considered. 

Other calculations have considered the multiple scattering in more detail, as diffusion of 
the electron. Migdal's 1956 calculation, in Sec. Ill, was the first; it has become something 
of a standard. A number of recent calculations of suppression are described in Sees. V to 
VII. 

Experimental work is surveyed in Sec. IX, beginning with the first cosmic ray studies 
of Landau-Pomeranchuk-Migdal (LPM) suppression shortly after Migdal's paper appeared. 
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These suffered from very limited statistics, and hence had limited significance. In the 1970's 
and 1980 's, a few accelerator experiments provided some data, but still with limited accuracy. 
In 1993, SLAC experiment E-146 made detailed measurements of suppression due to multiple 
scattering and photon-medium interactions. E-146 confirmed Migdal's formula to good 
accuracy, and at least partly inspired the more recent calculations in Sees. V to VII; because 
of this timing, some reference will be made to E-146 before the experiment is described in 
detail. 

Some more specialized topics will also be considered. Section X considers suppression in 
plasmas, where particles with similar energies interact with each other. Section XI shows how 
suppression mechanisms can affect electromagnetic showers, with a focus on cosmic ray air 
showers. Section XII surveys the application of LPM formalism to hadrons scattering inside 
nuclei, where color charge replaces electric charge, and gluons replace photons. Finally, a 
number of open questions will be presented. 

Earlier reviews (Feinberg and Pomeranchuk, 1956) (Akhiezer and Shul'ga, 1987) and 
monographs (Ter-Mikaelian, 1972) (Akhiezer and Shul'ga, 1996) have discussed these topics. 
A recent article (Feinberg, 1994) discusses the subject's early development. 

Unless otherwise indicated, energies will be given in electron volts (eV) and momentum 
in eV/c; He is 1.97 x 10~ 7 eV m and the fine structure constant a = e 2 /Tic ~ 1/137, where 
e is the electric charge. For both bremsstrahlung and pair conversion, E will represent 
electron/positron energy and k photon energy. 



II. CLASSICAL AND SEMI-CLASSICAL FORMULATIONS 



Landau and Pomeranchuk (1953a) used classical electromagnetism to demonstrate 
bremsstrahlung suppression due to multiple scattering; the suppression comes from the in- 
terference between photons emitted by different elements of electron pathlength. For those 
with quantum tastes, subsection II. B gives a simple semi-classical derivation based on the 
uncertainty principle. 



A. Classical Bremsstrahlung - Landau and Pomeranchuk 



The classical intensity for radiation from an accelerated charge in the influence of a 
nucleus with charge Z is 

d 2 I Z 2 e 2 u 2 f 

/ n x drexp {%[(ut - k • r(t)J) 



dudVt 



Att 2 c 3 



(6) 



where dr = v(t)dt, n is the photon direction and v(t) the electron velocity. The angular 
integration dVt is over all possible photon and outgoing electron directions. The classical 
bremsstrahlung spectrum can be found by assuming that v changes abruptly during the 
collision. Then, the integral splits into two pieces and the emission is easily found: 



d 2 I 


Z 2 e 2 


k x vi 


duidVl 


Atc 2 c 


k • vi — k c 



k x v 2 



k • v 2 — Iklc 



(7) 



where Vi and v 2 are the electron velocity before and after the interaction respectively. With 
small angle approximations, kv = u>(l — 1/27 2 )(1 — 9 2 /2), where # 7 is the angle between the 
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photon and incident electron direction and |k x v| = |k||v|# 7 , the bremsstrahlung angular 
distribution may be derived (for 7 >> 1, k <C E): 



dl 2 Z 2 e 2 7 4 |Av| 2 (l + 7 4 ^ 7 ) 



dudil tt 2 c 3 (1 + 7 2 # 2 ) z 



(8) 



where Av = vi — V2. The overall emission intensity can be found by integrating over dfl, 
and expressed in terms of the change in electron momentum q = 7ml Av|. Then, 

dl = 2Z 2 e 2 q 2 _ 
doj 37cm 2 c 3 

The complete bremsstrahlung cross section can be found by multiplying this by the elastic 
scattering cross section, as a function of q and converting from classical field intensity to 
cross section (Jackson, 1975, pg. 709): 

da 16Z 2 ar 2 fi™* dq 



dk 3k Jqmin Q 

where r e = e 2 /mc 2 = 2.8 x 1CT 15 m is the classical electron radius. The integral evaluates 
to In (q ma x / Qmin) , often called the form factor. This form factor is roughly equivalent to 
In (9 max /9 m i n ), where 9 is the electron scattering angle; this representation is convenient, 
because the distributions of angles depend on the suppression. The use of a, and implicitly 
h is needed to convert from field intensity to number of photons. 

The minimum and maximum momentum transfers correspond to maximum and min- 
imum effective impact parameters respectively. For a bare nucleus, the minimum impact 
parameter is the electron Compton wavelength A e = Ti/mc = r e /a = 3.8 x 1CT 13 m, so 
Qmax = 2mc. When the momentum transfer q max is larger than 2mc, the electron scatters 
by an angle larger than I/7. Then, coherence between the incoming and outgoing electron 
pathlength (implied in Eq. (7)) is lost and the cross section drops. The probability of such 
a large scatter is very small. This loss of coherence foreshadows how multiple scattering 
affects the electron trajectory and hence the radiation. The minimum momentum transfer 
occurs when qj_ = and <2W„ = q\\. 

However, nuclei are generally surrounded by electrons which screen the nuclear charge. 
For ultrarelativistic electrons, screening is important over almost the entire range of k. 
Screening shields the radiating electron from the nucleus at large distances, so the max- 
imum impact parameter is the Thomas-Fermi atomic radius a = 0.8Z _1 / 3 ao, for Z ^> 1 
(Tsai, 1974), where the Bohr radius a = h/amc = r e /a 2 = 5.3 x 1CT 9 m. Then, 
Qmin = 1-2571/ 'amcZ 1 / 3 and In (q max / ' q m in) = In (2.5/a;Z _1 / 3 ). A more detailed calculation 
finds ln(184Z- 1 / 3 ) for the form factor (Bethe and Heitler, 1934). 

This classical calculation is valid only for y — kj E <^ 1. A Born approximation quantum 
mechanical calculation gives a similar result, but covers the entire range of y = k/E (Bethe 
and Heitler, 1934): 

% - iJ w{»* + ^ + - < 184 ^ 1/3 » ^(* 2 + ^ + - *> 2 <»> 

where the radiation length 
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X = [Anar 2 e Z 2 In (I8AZ- 1 / 3 )]- 1 , 



(12) 



with n the number of atoms per unit volume. More sophisticated calculations, discussed in 
Section IV, include additional terms in the cross section and radiation length. With a more 
detailed screening calculation, a small constant may be subtracted from the form factor. 
However, Eq. (11) is a good semi-classical benchmark. 

If the interaction occurs in a dense medium, however, this treatment may be inadequate. 
The interaction is actually spread over the distance lf ] if the electron multiple scatters 
during this time, it can affect the emission. The multiple scattering changes the electron 
path, so that 

V (t) = v z (t) + v ± (t) = v(o)(i - e MS (t) 2 /2) + \v(o)\e MS (t)e (13) 

where is a unit vector perpendicular to the initial direction of motion and 0Ms{t) is 
the electron multiple scattering in the time interval to t. The particle trajectory is then 
r(t) = J v(t)dt, so that 

kT(t) = |v(0)| |k|(l-S) r ( i_^k )( ft + | k |0 f e MS {t)dt (14) 

2 Jo 2 Jo 

where # 7 is the angle between the photon and the initial electron direction. Landau and 
Pomeranchuk took 6 2 MS = (0 2 MS ). Over a distance d = v z t, the rms multiple scattering is 
(Rossi, 1952) 

(Ols) = (f) 2 ^ (15) 



where X is the radiation length and E s = mc 2 \J 'Att / a = 21.2MeV. 

The effect of multiple scattering can be found by inserting the trajectory of Eq. (13) 
into Eq. (6). This calculation simplifies with a clever choice of coordinate system. If the 
origin is centered on the formation zone, with vj_ = at z = 0, then the multiple scattering 
distributes evenly before and after the 'interaction'. Then, vx and v 2 are equally affected 
by multiple scattering, with 

k • Vl = k • v 2 = |k|(l - 1/2 7 2 )(1 - 6 2 /2)(l - 6 2 MS/2 /2) (16) 

and 



|k x Vl | = |k x v 2 | = |k|(l - l/2 7 2 )^ + ^ 5/2? (17) 

where 9ms/2 is the multiple scattering in half of the formation length. The angular addition 
in Eq. (17) is in quadrature because the angles are randomly oriented in the plane perpen- 
dicular to the electron direction. If the multiple scattering occurs on the same time scale as 
the interaction, then Eq. (8) becomes 



di 2 z 2 e 2 1 ^\ 2 \i + i\e 2 + e 2 MS/2 ) 2 } 



dudQ tt 2 c 3 [l+l 2 {0 2 + 2 MS/2 ) 



(18) 
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For 6 2 MS , 2 > I/7 2 + 2 , the radiation will be reduced and the angular distribution changed. 
If the multiple scattering in the formation zone is large enough, emission is suppressed. For 
9 7 <C I/7, this happens when 

Z?2 77i2 

^ < 44^- (19) 

Photons with lower energies will be suppressed. Equation (18) also shows that the an- 
gular distributions will be affected. Because multiple scattering can reduce the coherence 
length, a more detailed calculation is required to find the degree of suppression. Landau 
and Pomeranchuk (1953b) started with the expression 

dl Z 2 e 2 uo 2 r™ 



duj 4"7T 2 C 3 

which simplifies to 



/oo /*oo r 

(n x dri) ■ / (nx dr 2 ) dnexp{iio[(t 1 -t 2 ) -n-(r 1 -r 2 )/c]}, (20) 
-00 J— 00 J 



di z»a, r^ML -hiin + jy (21) 



where 



rfw 7rc 2 J-00 J-00 |r 12 



Ji = vi • v 2 - (Vl - ri2) 2 (V2 - ri2) f sin (g) + Sgcos^) -^) - 



and 



J 2 = (-2v 1 .v 2 )( gC ° sb) 9 7 Sinto) ). (23) 

Here g = w|r 12 |/c and r 12 = r(£i) — rfe)- The J\ term was evaluated by Landau and 
Pomeranchuk (1953b) while the J 2 term, neglected by Landau and Pomeranchuk was evalu- 
ated by Blankenbecler and Drell (1996), who found J\ = J 2 . Landau and Pomeranchuk still 
found the 'right' result, compensating by using a perpendicular momentum transfer due to 
multiple scattering twice the correct one. 

These integrals can be evaluated by using the electron trajectory given in Eqs. (13) and 
(14). Landau and Pomeranchuk gave a formula for the total (summed over time) emission. 



dl Z 2 e 2 uj r°° f°° dt ,,,, 
dT / — e lu}t 



d(J 7TC 2 J-oo J-oo t 3 

sin 



(* e MS (t)dt) 2 - te MS J* 0Ms(t)dt 



(24) 



where the integral over T = t\ — t 2 , will be factored out to get the emission per unit time. 
The remaining integral is evaluated by replacing the angular quantities by their average 
values. For example, 
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Landau and Pomeranchuk expressed their results in emission per unit time. In terms of the 
more prevalent total emission, 

dl AZ 2 e 2 r°° lv . , E 2 E 2 S X 2 . 
For uj <C E 2 c/E 2 X , the emission is 




dl 2itZ 2 e 2 m 2 c 3 uX 



du V 3 E,E 



(27) 



In the strong suppression limit, the field intensity dl/du ~ y/u, compared to an isolated 
interaction, where dl/duj is independent of u; the corresponding cross sections scale as 
1/y/uj and 1/uj respectively. The following subsection will describe a simple semi-classical 
derivation of the same result, and also discuss some of its implications. 



B. Bremsstrahlung - Quantum approach 

Bremsstrahlung suppression due to multiple scattering can also be found by starting with 
Eq. (5). Multiple scattering can suppress bremsstrahlung if it contributes significantly to q\\ 
(Feinberg and Pomeranchuk, 1956). Multiple scattering affects q» by reducing the electron 
longitudinal velocity. For a rough calculation, the multiple scattering angle can be taken as 
the average scattering angle. As before, the multiple scattering is divided into two regions: 
before and after the interaction. The longitudinal momentum transfer from the nucleus is 



911 = ^J(Eco S e MS/2 /c) 2 - (mc) 2 - \J[(E - k) cos6 MS/2 /c] 2 - (mc) 2 - k/c (28) 

where 9ms/2 is the multiple scattering in half the formation length, (E s /E) \Jlf/2X . The 
post-interaction scattering is based on the outgoing electron energy E — k. The inclusion 
of electron energy loss is a significant advantage of the quantum formulation. With some 
small angle approximations, 

km 2 c 3 k6 2 MS/2 
9 " = 2E(E - k) + (29) 

Multiple scattering is significant if the second term is larger than the first. This happens if 
Oms/2 > 1/7, or for 

k < k LPM = E{ *; ~ k) (30) 

&LPM 

where E LPM is a material dependent constant, given by 

^ m 4 c 7 X m 2 c 3 X Q a ,^ . . 

E LPM = = - 7.7 TeV/cm • X . (31) 

Table I gives E PP m in a variety of materials. This definition of E PP m was used by Landau 
and Pomeranchuk (1953a), Galitsky and Gurevich (1964), Klein (1993) and others, but is 
twice that used in some papers (Anthony, 1995, 1997) (Baier, 1996), and 1/8 of that in 
others (Stanev, 1982). The last choice is convenient for working with Migdal's equations 
(Sec. Ill), but not for the semi-classical derivation. 
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TABLE I. Suppression thresholds in a variety of materials. Listed are Z, Xq, density, -Elpm, Udie 
and E p . The carbon numbers is an average for graphite. BP stands for boiling point at atmospheric 



pressure. STP (standard temperature and pressure) is 20° C at 1 atmosphere. Standard rock, used 



to compare underground 


cosmic ray experiments, is defined as Z = 11, A = 


22 an 


dp = 


2.65g/cm 2 ; 


Xq is computed following 


f Tsaj 


H 974^1 Other 


mi in hpr^ a t*p 

J_l LLJ.11 UC1 □ CL1 C 


tttitti i Ti a rn pt t pi 

11 LJ111 1 1 J CL1 11C L L L I 


al, 1996). 




Material 


Z 


p (g/cm 3 ) 


Xq (cm) 


Elpm (TeV) 




Udie 


Ep 


hydrogen(liquid at BP) 


1 


0.07 


865 


6.6 PeV 


1.4 x 


10" 5 


4.6 PeV 


helium(liquid at BP) 


2 


0.125 


755 


5.8 PeV 


1.5 x 


10" 5 


4.3 PeV 


carbon 


6 


2.2 


19.6 


151 


5.5 x 


10" 5 


410 TeV 


aluminum 


13 


2.70 


8.9 


68 


6.0 x 


10" 5 


205 TeV 


iron 


26 


7.87 


1.76 


13.6 


1.0 x 


10" 4 


67 TeV 


lead 


82 


11.35 


0.56 


4.3 


1.1 X 


10" 4 


24 TeV 


tungsten 


74 


19.3 


0.35 


2.7 


1.5 x 


10" 4 


20 TeV 


uranium 


92 


18.95 


0.35 


2.7 


1.4 X 


10" 4 


19 TeV 


gold 


79 


19.32 


0.33 


2.5 


1.5 x 


10" 4 


19 TeV 


air(STP) 




0.0012 


30,400 


234 PeV 


1.3 x 


10" 6 


15 PeV 


water 




1.00 


36.1 


278 


3.9 x 


10" 5 


540 TeV 


std. rock 


11 


2.65 


10.0 


77 


6.0 x 


10" 5 


230 TeV 



For k E, Eq. (30) matches the classical result: photon emission for k < k^pu = 
E 2 /Elpm is reduced. For example, for 25 GeV electrons in lead, photon emission below 
125 MeV is suppressed. For higher energy electrons, a large portion of the spectrum can be 
suppressed. For 4.3 TeV electrons in lead, fewer photons with k < 2.15 TeV are radiated. 
At the endpoint of the spectrum, k — > E, lf approaches zero, so the Bethe-Heitler cross 
section always applies there. 

Different calculations have found slightly different coefficients in Eq. (30). Landau and 
Pomeranchuk match Eq. (30), but Feinberg & Pomeranchuk (1956) and the introductions 
to Anthony (1995, 1997) used the simpler criteria (9ms) > 1/7, and find a suppression twice 
that given here. Use of the minimum uncertainty principle ApAx > K/2 (Schiff, 1968, pg. 
61) in Eq. (5) would shorten Ifo and reduce the suppression. 

Because 9ms/2 both depends on If, and also partly determines If, finding the suppression 
requires solving a quadratic equation for If. 
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If multiple scattering is small, this reduces to Eq. (5). Where multiple scattering dominates 



The bremsstrahlung cross section scales linearly with the distance over which coherence 
is maintained, or the formation length. It is convenient to define a suppression factor S, 
giving the emission relative to Bethe Heitler: 



da/dk lj_ \ kE LPM ( * 

da BH /dk l f0 \E(E-k) 1 } 

and the dN/dk ~ 1/k found by Bethe and Heitler (1934) changes to dN/dk ~ 1/y/k. 

Figure 1 compares bremsstrahlung cross sections for 25 GeV and 1 TeV electrons in 
lead. The Bethe-Heitler cross section is compared with three approaches to suppression: 
the full semi-classical suppression Eq. (32), the prevalent low energy limit Eq. (34), and 
Migdal's calculation, included as a standard. Equation (32) predicts considerably more 
suppression than Migdal. Numerically, Eq. (34) is closer to Migdal's results, but the required 
approximation is unjustified in the transition region k ~ E[E — k)/ Elpm- Better agreement 
would be found with a larger Elpm, as would be given by the minimum uncertainty principle. 
With only S considered, it is impossible to separately determine Elpm and the overall cross 
section normalization. The onset of suppression must be seen to separate these two factors. 

Nevertheless, Eq. (34) demonstrates some implications of bremsstrahlung suppression 
(Landau and Pomeranchuk, 1953b). With suppression, the number of photons emitted 



per radiation length is finite, scaling as ^E/E L pm for E > E L pm- Electron dE/dx due 
to bremsstrahlung is also reduced; instead of rising linearly with E, it is proportional to 
\JEE L pm- Figure 2 shows the relative bremsstrahlung energy loss, [dE j dx)LPM for Migdal 
compared with the Bethe-Heitler prediction. 

Besides the photon spectrum, LPM suppression also affects the photon angular distribu- 
tion. With the photon emission angle 6 1 included, 

9|| — Pe — p' e ~ k/ccos9 y = \J (E / c) 2 — (mc) 2 — \J[(E — k)/c} 2 — (mc) 2 — k/ccosO^. 

(35) 

The change in electron direction is assumed to be negligible. Then, 

m 2 c 3 k k6 2 
g '' ~ 2E(E - k) + ~c (36) 

For k E, If = Ifo/il+^O 2 ). This formula may be used to derive the angular distribution 
of bremsstrahlung photons (Landau and Pomeranchuk, 1953b). When 6 1 increases q\\, the 
multiple scattering term becomes less important, so there is less suppression for 7# 7 > 1. 
With both multiple scattering and a finite 9 1 , 
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(a) 



125 



250 
k (MeV) 



375 



500 




500 750 1000 

k (GeV) 

FIG. 1. Comparison of the differential energy weighted cross sections per radiation length, 
X^nkda/dk for Bethe-Heitler radiation (solid line), Migdal's detailed suppression calculation (dot- 
ted line), quadratic suppression, Eq. (32) (short dashes), and the strong suppression limit, Eq. 
(34) (long dashes), for (a) 25 GeV electrons and (b) 1 TeV electrons incident on a lead target. 
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0.8 



75 0.4 




10" 



10° 



1 10* 10* 

Electron or Photon Energy/E LPM 

FIG. 2. The relative electron energy loss dE/dx(Migd&l)/dE/dx(Bethe-lleitler) from 
bremsstrahlung (solid line) and the relative photon conversion cross sections d(Migdal)/cr(Bethe- 
Heitler). Both curves are for lead, with electron/photon energy in units of Elpm = 4.3 TeV. 



2hE(E - k) 
km 2 c 3 



1 + 



E 2 l 



f 



2m 2 c 4 X 



+ 



6 2 y E(E-k)T-> 



m 2 c 4 



(37) 



Suppression is large when the multiple scattering term is larger than the other terms. Then 



S(9) 



kE LPM (l + 7 2 9 2 ) 
\ E(E - k) 



(38) 



Suppression disappears rapidly as 9 1 rises. 

When 5(0) <C 1, the angular distribution is broadened. It follows from Eq. (38) that 
multiple scattering broadens the angular distribution from (0 7 ) ~ I/7 to (9 y ) ~ \/ r y\/~S 
(Landau and Pomeranchuk, 1953b). Starting from a different tack, Galitsky and Gurevich 
(1964) found that (0 7 ) is determined by the electron multiple scattering over the distance 
//, giving the same algebraic result. 

This increase in (9 y ) is difficult to observe. The angular distribution of bremsstrahlung 
photons in a thick target is dominated by changes in electron direction due to multiple 
scattering before the bremsstrahlung occurs. However, for sufficiently thin targets, multiple 
scattering will be small, so that, if S is small enough, then the photon emission angles 
dominate over multiple scattering. For thin targets, the photon spectrum measured at 
angles 9 1 3> I/7 should exhibit less suppression than at smaller angles. 
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C. Photon Interactions with the Medium 



Ter-Mikaelian (1953a,b) pointed out that photon interactions can also induce suppres- 
sion. Photons can interact with the medium by coherent forward Compton scattering off 
the target electrons, producing a phase shift in the photon wave function. If this phase shift 
is large enough, it can cause destructive interference, reducing the emission amplitude. Ter- 
Mikaelian used classical electromagnetism in his analysis, calculating suppression in terms 
of the dielectric constant of the medium, 

e(k) = 1 - {huj p ) 2 /k 2 (39) 



where tu p , the plasma frequency of the medium, is y AnnZe 2 /to. This is equivalent to giving 
the photon an effective mass hu p /c 2 . The relationship between k and p 7 becomes p 7 c = y/ek, 
so 

«l=W.-W?/c=^ + ^. (40) 



The formation length is then 



U = (41) 



k 2 + k 2 p 



where k p = •yhuip. When dielectric suppression is strong, q\\ is dominated by the photon 
interaction term and If becomes independent of E: If = 2ck/fiui 2 . As with LPM suppression, 
the cross section is proportional to the path length that can contribute coherently to the 
emission, so S is the ratio of the in-material to vacuum formation lengths: 

h 2 

(42) 



k 2 + k 2 p 



For k < k p , bremsstrahlung is significantly reduced. This happens for y = k/E < y die , where 

y die = hujp/mc 2 (43) 

is a material dependent constant. For lead, hw p = 60 eV, so y die ~ 10~ 4 . Table I lists y die 
for a variety of materials. 

The same result can be obtained classically by including the dielectric constant in Eq. (1), 
so $ = exp {ikt[l — \v\y/ecos(9 1 )/c]/h}. The case cos (# 7 ) = c/\v\y/e, which gives an infinite 
formation length, corresponds to Cerenkov radiation (Ter-Mikaelian, 1972, pg. 196). 

Because photon emission angles are determined by the kinematics, a finite # 7 affects 
dielectric suppression the same way as it does LPM suppression. Including 6> 7 , 

= 2ftc/c 7 2 

f P(l + 7 2 ^) + A;2- [ > 

For large suppression and 7# 7 > 1, S — (kO^/hujp) 2 is independent of E. In this limit, the 
angular spread is < # 7 >~ hu p /k (Galitsky and Gurevich, 1964). 
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Ter Mikaelian (1972, pg. 127) pointed out that the dielectric constant also affects q, 
in the form factor logarithm. The complete cross section for k <C E is then 



da _ lQZ 2 ar 2 e k 
dk ~ 3k 2 , 



In (l84Z- 1 /3^i + (^)2 ) 



(45) 



Because k 2 ~ n, except for the logarithmic term, photon emission is independent of the 
density! As the density rises, increasing the number of scatters, suppression rises in tandem, 
leaving the total photon production constant. 

This suppression is sometimes known as the longitudinal density effect, by analogy with 
the transverse density effect (Jackson, 1972, pg. 632) which reduces ionization dE/dx. It 
is also known as dielectric suppression. Unfortunately, a quantum mechanical calculation 
of dielectric suppression has yet to appear, nor has dielectric suppression been described in 
terms of Compton scattering. 

Because dielectric suppression and the LPM effect both reduce the formation length, 
the effects do not merely add; the total q\\ must be calculated, and from that If and the 
suppression can be found. Feinberg and Pomeranchuk (1956) showed that when k p > k^pM 
(i.e. E < UdieELPM), then dielectric suppression overwhelms LPM suppression, and only the 
former is observable. For higher electron energies, LPM suppression is visible for 



D. Bremsstrahlung Suppression due to Pair Creation 

Landau and Pomeranchuk (1953a) pointed out that, at the highest energies, If can 
approach a radiation length. Then, the partially created photon can pair create part way 
through the formation zone. This destroys the coherence between different parts of the 
formation zone, reducing the amplitude for photon emission. Unfortunately, there has been 
little attention to this problem, and the available results are quite crude. 

For k ^> mc 2 , the pair creation cross section is independent of energy: a 1 = 
(28/9)ar2Z 2 ln(184Z~ 1 / 3 ). This constant cross section limits the formation length to 
roughly X . Neglecting other suppression mechanisms, Ifo > X when k < 2HE(E — 
k)/X m 2 c 3 . However, dielectric suppression and LPM suppression limit the range of ap- 
plicability. With both these mechanisms considered, pair creation further reduces photon 
emission when (Galitsky and Gurevitch, 1964): 



The coefficients given here differ slightly from the original because Galitsky and Gurevitch 
used a slightly different approach from that presented here. With other factors considered, 
this mechanism is visible for 



k > k cr = (k 4 p /k LPM ) 1/3 . 



(46) 




(47) 



E> E p = 



V2c 



(48) 
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FIG. 3. Schematic view of bremsstrahlung da/dk when several suppression mechanisms are present. 
For E < E p , the pair creation suppression disappears and LPM suppression connects with dielectric 
suppression. 



is a material dependent constant. For E ~ E p this mechanism dominates in a narrow 
window around k = Tioj^Xq/c] for E ^> E p , the range is k p _ = hu p X /2c < k < k p+ = 
AhcE(E — k)/(X Eg), as is shown in Fig. 3. In this region, the suppression factor is 



Xo =J n^kX^_ 
l f0 'Hi I ■'.(!'. - k) 



and da/dk is independent of k. For k < k p ^, dielectric suppression dominates, while for 
k > k p+ , LPM suppression is dominant. E p ranges from 25 TeV for lead to 15 PeV for sea 
level air; other values are given in Table I. For lead, when E ^> E p , the photon 'window' is 
1.2 xlO 7 eV< k < 1.6xl(T 19 £ 2 (eV), while for air it is8.3xl0 7 eV < k < 3.0 x l(T 24 £ 2 (eV). 

Similarly, bremsstrahlung of a sufficiently high energy photon can suppress pair produc- 
tion. The bremsstrahlung can affect the overall pair production rate if the emitted photon 
contributes significantly to q\\ of the entire reaction. 

These formulae are only rough approximations. At high enough energies, the pair cre- 
ation cross section is itself significantly reduced because of LPM suppression and X in 
Eqs. (47)- (49) should be increased to account for this. As Fig. 2 shows, for k 3> E LPM , the 
pair conversion length rises significantly. Then, bremsstrahlung and pair creation suppress 
each other and the separation of showers into independent bremsstrahlung and pair creation 
interactions becomes problematic. For this, a new, unified approach is needed. 



E. Surface effects and Transition Radiation 



The discussion far has only considered infinitely thick targets. With finite thickness 
targets, the effects of the entry and exit surfaces must be considered. At first sight, this 



16 




FIG. 4. Diagram of transition radiation caused by multiple scattering. The pancaked electromag- 
netic fields are broadened by multiple scattering from their free space width 1/7 to + MS ^ 2 , 
where 6ms depends on If, and hence on k. 



appears straightforward: the only effect being a reduction in the multiple scattering when 
the formation zone sticks out of the target, and hence there is less suppression. However, 
in addition to reduced suppression, multiple scattering produces a new kind of transition 
radiation. 

Conventional transition radiation occurs when an electron enters a target, and the elec- 
tromagnetic fields of the electron redistribute themselves to account for the dielectric of the 
medium. In the course of this rearrangement, part of the EM field may break away, be- 
coming a real photon. Because transition radiation has been extensively reviewed elsewhere 
(Artru, Yodh and Mennessier, 1975), (Cherry, 1978), (Jackson, 1975), it will not be further 
discussed. However, the formula for radiation, neglecting interference between nearby edges, 
is given here for future use. The emission is (Jackson, 1975, pg. 691): 



dN _ a 
dk nk 



2k 2 \ . hi. 



(50) 



photons per edge. 

The additional transition radiation occurs because multiple scattering changes the trajec- 
tory of the electron. The variation in electron direction widens the directional distribution of 
the electromagnetic fields carried by the electron, as is shown in Fig. 4. Scattering broadens 
the EM fields from their free space width I/7 to a If (and hence k) dependent value. As the 
EM field enters the target and realigns itself, it can emit transition radiation. 

Classically, this new form of transition radiation is closely related to the old. Both 
depend on the difference in If in the two materials (Ter-Mikaelian 1972, pg. 233), with the 
complete radiation 

dN - 20W VW (51) 



dkdO^ -kTl 6 



where If and l'f are the formation lengths in the two media. If If = If, then there is no 
transition radiation. Conventional transition radiation can be derived from this formula 
by focusing on the dielectric constant of the media, while multiple scattering transition 
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radiation can be calculated by focusing on 6ms- Of course, the complete spectrum includes 
both contributions. 

The energy spectrum is given by integrating over # 7 ; the integral is complicated because 
the maximum # 7 depends on k. For E > E LPM and k > k p , (Ter-Mikaelian, 1972, pg. 235) 



dN «/ 1 + Jl + Ak LPM /k 2 \ 

~ — In ( ) H 1 ) 

dk ~k V 2 1 • \ 1 • 1/.7./..V/7.- 7 



(52) 



per edge. The total radiation may be found by integrating dN/dk up to &lpm- Since 
khPM ~ E 2 , the total energy lost by the electron AE rises as E 2 , in contrast to conventional 
transition radiation, where the loss is proportional to E. When AE becomes a significant 
fraction of E, quantum effects must become important. 

These calculations assume that the transition is instantaneous, neglecting the sharpness 
of the surface. They also neglect coherence between nearby edges. More quantitative es- 
timates are discussed in Sec. IV. C. For pair creation, one expects similar surface effects; 
unfortunately these have yet to be worked out. 



F. Thin Targets 

In extremely thin targets, neither dielectric effects nor multiple scattering produce enough 
of a phase shift to cause suppression. Suppression due to multiple scattering disappears 
when the total scattering angle in the target is less than 1/7. This occurs for targets with 
thickness T < (mc 2 /E s ) 2 X pa X /1720. With dielectric suppression, If — > as k — > 0, so 
some dielectric suppression remains for any T, for photon energies k < Thu 2 /2c. Of course, 
for extremely thin targets transition radiation will dominate over bremsstrahlung. 

For multiple scattering, intermediate thickness targets with X /1720 < T < lf are of 
intererest because the entire target acts as a single radiator. The emission can be found 
from the probability distribution for the total scattering angle in the entire target, either 
classically (Shul'ga and Fomin, 1978) or with quantum mechanical calculations (Ternovskii, 
1960). These calculations are more complex than those presented earlier, because they use 
a distribution of scattering angles rather than a single average scattering angle. The total 
radiation from the target is (Ternovskii, 1960) 

where ( = 7#/2, and 9 is the scattering angle. Since 9 is independent of k, this formula 
has the same k dependence for k <C E as the Bethe-Heitler calculation. The scattering can 
be treated as Gaussian distribution, with rms scattering angle 9 , where, for thin targets 
(Barnett et ai, 1996) 

o = §J|^(l + O.O381n-p. (54) 



For the relevant range of thicknesses, neglecting the logarithmic term changes the suppression 
factor by at most a few percent. 
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FIG. 5. Suppression factor Nt/Tcjbh for thin targets. The solid line follows Eq. (53) for E=25 
GeV. The factor is only slightly energy dependent, as long as T > Xo/1720 and T < If. The 
dashed curve is the logarithmic approximation, Eq. (55). 

Eq. (53) can be evaluated numerically. For very thin targets, it matches the Bethe- 
Heitler spectrum, except for a factor of 4/3. For thicker targets, the suppression factor, 
S = Nt/Tctbh is shown in Fig. 5. 

In the limit T > (mc 2 / E S ) 2 X , (but with T < l f0 ), dN/dk ~ InT (Ternovskii, 1960). 
The intensity varies logarithmically with the target thickness! A slightly more detailed 
calculation finds (Shul'ga and Fomin, 1996) 



This approximation, shown by the dashed line in Fig. 5, overestimates Eq. (53) by 10-20%. 

Shul'ga and Fomin (1998) found that the use of a screened Coulomb potential, rather than 
a Gaussian scattering distribution, leads to an equation like Eq. (55), but with additional 
terms. This calculation shows the effects of scattering on an electron by electron basis. The 
more an electron scatters in a target, the higher the average radiation. 

A similar expression should apply for pair creation - the formation length is the same. 

G. Suppression of Pair Creation 

Multiple scattering can also reduce the cross section for 7 — > e + e~. The relationship 
between pair creation and bremsstrahlung, Fig. 6, is clear, and the two Feynman diagrams 
easily map into each other. The crossing does change the kinematics of the process. Since 
it is the electron and positron that multiple scatter, and they must have energies lower than 




(55) 
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that of the initial photon, suppression occurs only at higher incident particle energies. To 
avoid confusion, k will continue to refer to photon energy, with the produced pair having 
energies E and k — E. 

Pair production is not possible classically. However, Landau and Pomeranchuk (1953b) 
give some simple arguments why pair creation should be sensitive to its environment. The 
momentum transfer for pair production is (Feinberg and Pomeranchuk, 1964) 

9 „ = k/c - y/(k - EY/c* - m*<? - y^/c* - mV = ^j^y • (56) 

Because q\\ is unchanged when E and k — E are interchanged, either can represent the 
electron or positron. The formation length can be expressed in terms of the two final state 
momenta or in terms of the invariant mass of the created pair. Not surprisingly, lf for pair 
production is similar to the bremsstrahlung case: 

//0 = m^k ■ (57) 

It might seem surprising that k is in the denominator of Eq. (57). But, lf becomes a 
maximum for E = k — E = k/2; then lf = hk/2m 2 c 3 , and lf rises with k. If E <C k, 
then this equation reduces to lf = 2hE/m 2 c 3 and Z/ is ver y short. This asymmetric energy 
division corresponds to a pair with a large invariant mass. In terms of pair mass, M p , 

One difference between pair creation and bremsstrahlung is that the multiple scattering now 
applies to two separate particles. The lower energy particle scatters more, and so dominates 
the additional qn. The scattering is taken over lf/2, as if the charged particles are produced 
in the middle of the formation zone, and the result is very similar to Eq. (32): 



If = I 



fo 



1 + 



EH 



-i 



2m 2 c 4 X . 

where Ifo now refers to Eq. (57). When the second term is dominant, 



(59) 
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/ kE LPM 2hk IE L pm , rn , 

similar to Eq. (32). Then, 

a / E LPM k M p E LPM 

S= iE(k—Ej = ^ — - (61) 

For a given fc, the suppression is largest when E m (k — E) m k/2. There is no suppression 

for E « or E pa fc. For E 3> E LPM , the total cross section scales as ^ E LPM /k. Figures 2 
and 10 show how S 1 drops as fc rises. 

As with bremsstrahlung, the emission angles can affect suppression. The relevant angular 
variables are 9 e + and 9 e -, the angles between the outgoing particles trajectories and the 
incoming photon path. If either angle is larger than k/mc 2 , then the formation length is 
shortened and suppression reduced. 

Because of the high photon energy, there is no apparent analogy to dielectric suppression 
for pair creation. 



H. Muons and Direct Pair Production 

Electromagnetic processes involving muons can also be suppressed. However, because 
the muon mass m M ^> m e , the effects are much smaller. For a fixed energy, the formation 
length is reduced by (m e /m M ) 2 ~ 1/40,000. For muons, Eqs. (30) and (34) hold, but with 
Elpm replaced by 



Elpmm = m £ E f° « 1-38 x 10 22 eV/cm • X . (62) 



a 



This energy is high enough that LPM suppression is generally negligible for muon 
bremsstrahlung and pair creation. 

For muons, dielectric suppression still occurs for y < hup/m^c 2 , about 10 -7 in solids. 
This is very small, but perhaps not unmeasurable. 

Unlike electrons, high energy muons have a significant cross section for direct pair pro- 
duction, fi~N — > ji~e + e~N . This process is similar to bremsstrahlung followed by pair 
creation, except that the intermediate photon is virtual. Both the \i and the final state elec- 
trons are subject to multiple scattering. The formation length can be calculated by treating 
the pair as a massive photon, starting from 

m 2 c 3 k M 2 c 3 



g|| = ^ E2 ~ m V ~ ^ E - k f - W - V^ 2 - M p c " = 2E (E -k) + it (63) 

where here E is the muon energy, k is the virtual photon energy and M p is the pair 
mass. Compared to bremsstrahlung, If is only decreased for k/E < M p /m^. The final 
state has three particles that can multiple scatter. Since the incoming muon is very en- 
ergetic, it exhibits little multiple scattering. The electron and positron are less energetic, 
and multiple scatter more; the contribution to q\\ due to (both) their multiple scattering is 
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E 2 s M 2 lf/Am 2 kX . This is significant (the cross section is reduced) for If > (E s /mc 4 ) 2 X /2 
when k/E < M p /m p . Suppression is easiest for symmetric pairs because they lead to the 
longest If, but, even then, energies above 10 17 eV are required to observe suppression. 

Although it is much less probable, electrons can also lose energy by direct pair production, 
e~N — > e + e~e~N. For electrons, M p > m, so the second term in Eq. (63) always dominates, 
and // = 2hk/M 2 c 3 . As with muons, suppression occurs for If > (E s /mc 2 ) 2 X /2, albeit 
without the restriction on k/E. Still, suppression requires energies almost as high as the 
muon case. 

The lack of suppression for direct pair production is an initial demonstration that higher 
order diagrams typically involve larger q\\ than simpler reactions. Therefore, higher order 
processes are less sensitive to their environment, and, when suppression is large, higher order 
diagrams become more important. 



I. Magnetic Suppression 

External magnetic fields can also affect the electrons trajectory, and hence its radiation. 
This section will consider the effect of the change in electron trajectory on bremsstrahlung 
emission, neglecting the closely connected synchrotron radiation emitted by the same field. 

An electron will be bent by an angle 

Ap _ eBlf sm(j) B 

e B/ 2 - — - — (64) 

in a distance lf/2 in a uniform magnetic field B. Here, 0b is the angle between the elec- 
tron trajectory and the magnetic field. As with multiple scattering, if 6 B / 2 > 1/7, then 
bremsstrahlung is suppressed. This happens when (Klein, 1993) 

lBsm</) B 

y <Vb = ^ (65) 



where B c is the critical magnetic field, B c = m 2 c 3 /eh = 4.4 x 10 13 Gauss. 

The bending angle 9 B / 2 accumulates linearly with If, in contrast to the LPM case where 
&MS/2 ~ l 2 f] this leads to a stronger k dependence than with LPM scattering. If 9 B / 2 is 
treated in the same manner as MS / 2 in Eq. (28), then 



l f = I 



fo 



1 + 



/ mcB sin 4> B lf\ 2 



V HB, 



c 



(66) 



This is a quartic equation for If, compared with the quadratic found with multiple scattering. 
In the limit of strong magnetic suppression (If <^ lf ), the suppression factor /////o has a 
form similar to the LPM effect (Klein, 1997): 

where E B = mc 2 B c /B sin (p B . Figure. 7 shows the suppression for three different values of 
E/E B . 
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FIG. 7. Suppression factor 5 for magnetic suppression, for E = O.IEb (long dashes), E = IEb 
(solid line) and E = WEb (short dashes). 

Because magnetic suppression has a weaker k dependence than dielectric suppression, 
it is only visible when dielectric suppression does not apply, i.e. for E > ydieEs, and then 
for db > y > Vdie- For 25 GeV electrons in saturated iron (~ 20 kG), Eb ~ 15 PeV 
and db ~ 2.5 x 10~ 5 , so the magnetic effect will be hidden. At higher electron energies, it 
becomes quite visible. For a 1 TeV electron in a 4T field, as will be found at in the CMS 
detector at LHC, photons with energies below 900 MeV are suppressed. 

Suppression should also occur for pair production. A similar calculation finds 



For symmetric pairs (maximum suppression), S = (AEs/k) 2 ^ 3 . Because the magnetic bend- 
ing is quite deterministic, in contrast to multiple scattering which is statistical, this semi- 
classical calculation may be more accurate than that for multiple scattering. 

Baier, Katkov and Strakhovenko (1988) considered bremsstrahlung suppression in a mag- 
netic field, for both normal matter (screened Coulomb potentials) and e + e _ colliding beams. 
They used kinetic equations to find the radiation to power law accuracy, in both strong and 
weak field limits. These results are similar to Eq. (67). 

When the magnetic field is confined to the material, magnetic suppression should also 
produce transition radiation. This should be most visible with ferromagnetic materials. 
Equation (51) could be used to find the spectrum. 

Both the semi-classical calculation and the more accurate result neglect synchrotron 
radiation. Because the formation length scales for bremsstrahlung and synchrotron radiation 
are similar, they are important in the same kinematic regions. A complete calculation should 




(68) 
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treat them together, calculating the electron trajectory due to the combined field, and then 
calculating the radiation for that trajectory. 



J. Enhancement and Suppression in Crystals 

So far, we have considered only amorphous materials. In crystals, however, the regu- 
larly spaced atoms produce a huge range of phenomena, because the interactions with the 
different atoms can add coherently (Williams, 1935). When the addition is in phase, en- 
hanced bremsstrahlung or pair production results, while out of phase addition results in a 
suppression. 

Electrons can interact with atoms with impact parameters smaller than h/q±. The rela- 
tive phase depends on the spacing between atoms measured along the direction of electron 
motion. The phase difference for two interactions is $ = exp {i[ut — k • a(£)]} where a is 
the atomic position. If $ has the same phase for two nuclei, then the emission amplitudes 
add coherently. 

Including 7 , the phase difference between two adjacent sites is 



where a is the spacing between two atoms along the direction of electron motion. If the 
nuclei are spaced so that 



where n is an arbitrary integer, then $ = 1 and the addition is coherent. As with Cerenkov 
radiation, for certain # 7 , the phase is always zero, and If — > oo, implying infinite emission 
(from an infinite crystal). Conversely, if a/ccos(# 7 ) = (2n + l)2irhE(E — k)/m 2 c 3 , there is 
complete destructive interference. 

The large set of variables in Eq. (70) gives rise to a variety of effects. As the incident 
electron direction (affecting a) and/or # 7 vary, the interference will alternate between con- 
structive and destructive, producing peaks in the photon energy spectrum for most sets of 
conditions. 

Although Eq. (69) predicts infinite coherence, in a real crystal several factors limit the 
coherence length. One of these is the thermal motion of the atoms. When the rms thermal 
displacement of the atoms is larger than l/q\\, the coherence is lost. When Q 1 is large, the 
transverse separation between the electron and photon can limit the coherence. 

Changes in the electron trajectory can also reduce the coherence length. The crystalline 
structure can generate very high effective fields, causing strong bending, known as channeling 
(S0rensen, 1996); this bending can limit the coherence length. Multiple scattering can also 
change the electron direction, and limit the coherence (Bak et ai, 1988). Finally, crystal 
defects and dislocations can also limit coherence. Because this is a vast subject, with several 
good reviews available (Palazzi, 1968) (Akhiezer and Shul'ga, 1987) (Baier, Katkov and 
Strakhovenko, 1989) (S0rensen, 1996), this article will not consider regular lattices further. 



$ = exp — k ■ v]acos(# 7 )} 



(69) 




4imhE(E - k) 



(70) 
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K. Summary & Other suppression mechanisms 

The suppression mechanisms discussed so far are summarized in Table II, in order of 
increasing strength. Many other physical effects can lead to suppression of bremsstrahlung 
and pair production. Many of them involve partially produced photon interactions with the 
medium. Some of the interactions that can affect bremsstrahlung are photonuclear interac- 
tions, real Compton scattering, and, at lower energies, a host of atomic effects including K 
and L edge absorption and a variety of optical phenomena. 



TABLE II. Summary of different bremsstrahlung suppression mechanisms. 



Region 


Source 


dN/dk scaling 


Maximum k 


Bethe-Heitler 




k- 1 


E 


LPM 


multiple scattering 




kLPM 


Magnetic 


magnetic field 




VbE 


Pair Creation 


pair creation 


k° 


k p+ 


Dielectric 


photon interactions 


k 1 


VdieE 



Photonuclear interactions can have an effect similar to that of pair conversion - the partly 
created photon is destroyed. Because the photonuclear cross section is much smaller than 
the pair conversion cross section, this is a small correction. Real Compton scattering can also 
effectively destroy the photon. Toptygin (1964) treated these reactions as imaginary (ab- 
sorptive), higher order terms to the dielectric constant of the medium. The bremsstrahlung 
plus real Compton scattering of the partially produced photon is in some sense a new class of 
radiation, with its own Feynman diagram. Because Compton scattering involves momentum 
transfer from the medium, // is short, and so, in some regions of phase space (when dielec- 
tric suppression is large), Toptygin found that this diagram can be the dominant remaining 
source of emission. 

Other photon absorption mechanisms occur at lower photon energies. For example, K or 
L edge absorption produces a peak in the photon absorption spectrum. If, over a formation 
length, the absorption probability due to these peaks is significant, then suppression can 
occur. These peaks are difficult to observe because of competition from transition radiation, 
which is enhanced at the same energies (Bak et ai, 1986). For optical photons, a host 
of atomic effects can affect the dielectric constant. These variations can also introduce 
suppression (Pafomov, 1967). 

III. MIGDAL FORMULATION 

More quantitative calculations are more complex. The first quantum calculation, by 
Migdal (1956, 1957) is still considered a standard. Migdal treated the multiple scattering 
as diffusion, calculating the average radiation per collision, and allowing for interference 
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between the radiation from different collisions. When collisions occur too close together, 
destructive interference reduces the radiation. 

This approach replaces the average multiple scattering angle used earlier with a realistic 
distribution of scattering, and, hence, of of the electron path. Migdal also allows for the 
inclusion of quantum effects, such as electron spin and photon polarization. 

Migdal treated multiple scattering using the Fokker-Planck technique (Scott, 1963). This 
technique is used to solve the Bolzmann transport equation for F(0,z), where F is the 
probability distribution of particles moving at an angle 9 with respect to the initial direction. 
The equation is: 



dz 



= -lo (z)F(9,z) + 9'd9' d/3'W(x,z)F(e',z) (71) 

JO Jo 



where ujq is the scattering probability per unit thickness, given by the integral of W(9, z) 
over all angles. Here, 9 is the scattering angle and j3' is the azimuthal angle. W(x, z) is 
the single scatter angular distribution. The angle x is the angular opening between vector 
representing 9 and 9': x 2 = 9 2 + 9' 2 — 299' cos f3'. The W dependence on z allows for 
inhomogeneity in the material; otherwise W is independent of z. For each of the trajectories 
allowed by the diffusion, Migdal calculated the photon radiation, including electron spin and 
photon polarization effects. 

The Fokker Planck method is valid if W(x, z) is sufficiently sharply peaked at % = 0, so 
that it has a finite mean square (x 2 )z and that F(0', z) — F(9, z) can be accurately approxi- 
mated by a second order Taylor expansion in 9' — 9. Some calculations (Scott, 1963) lead to 
a Gaussian distribution for F, with a mean multiple scattering angle (x 2 )z- Unfortunately, 
a Gaussian distribution underestimates the number of scatters at angles larger than a few 
times 9q (Eq. (54)). This problem limits the accuracy of this calculation. The problem is 
somewhat exacerbated because If is relatively short, so the number of scatterings is fairly 
small. 



A. Bremsstrahlung 

With these calculations, updated with a more modern form factor, the Migdal cross 
section for bremsstrahlung is 

*%L = *2gW + 2[1 + (! - » m .))Z> In (|i) (72) 

where G(s) and <f>(s) are the suppressions of the electron spin flip and no spin flip portions 
of the cross section respectively, 

w V4 2 Jo sinh(t/2) J v 1 

0(s) = 12s 2 ( J™ e~ sx coth (x/2) sin (sx)dx^j - Qns 2 . (74) 

The factor £(s) is 
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FIG. 8. Migdals G(s) (dashed line) and <f>(s) (solid line). 4> governs the suppression for y <C 1. 



with Si 



£(s) = 2 

£( S ) = l + ln( S )/ln( Sl ) 

m = i 

Z 2 / 3 /184 2 . Here, 



< si) 

(si < s < 1) 

(S>1) 



Elpm^ 



8E(E-k)as)' 



(75) 



(76) 



Migdal gave infinite series solutions for G(s) and 0(s). However, they may be more simply 
represented by polynomials (Stanev et al, 1982): 

(j)( s ) = 1 - exp ( - 6s[l + (3 - tt)s] + s 3 /(0.623 + 0.796s + 0.658s 2 )) 

^(s) = 1 - exp ( - 4s - 8s 2 /(l + 3.96s + 4.97s 2 - 0.05s 3 + 7.5s 4 )) 

G(s) = 3^(s) -20(s). (77) 

These functions are plotted in Fig. 8. 

For k <C E, s ~ 1/(7#ms)- For s » 1, there is no suppression, while for s « 1, the 
suppression is large. In the absence of suppression G(s) = <f>(s) = 1, and Migdal matches 
the Bethe-Heitler cross section. For small suppression, where s is large, G(s) = 1 — 0.22/s 4 
and (f>(s) = 1 — 0.012/s 4 . For strong suppression, where s — > 0, G(s) = 12ns 2 and 0(s) = 6s. 

If a Coulomb scattering distribution is fit with a Gaussian, the mean of the Gaussian 
grows slightly faster than yJTf. This increase is reflected in the the logarithmic rise in £. For 
sufficiently high energies, q m ax is limited by the nuclear radius Ra, which Migdal approx- 
imated Ra = Q.har e Z 1 ^ . The nuclear radius limits q max to Ti/Ra, at higher momentum 
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transfers the nuclear structure is important and Eq. (72) loses accuracy. The fortuitously 
simple '2' coefficient comes from the chosen approximation for This cutoff is only 

reached for extremely large suppression, and is hence of limited importance; higher order 
terms are likely to dominate at this point. 

One difficulty with these formulae is that £ depends on s, which itself depends £, so 
the equations must be solved recursively. To avoid this problem, Stanev and collaborators 
(1982) developed simple, non-iterative formulae to find s and £(s). They removed £(s) from 
the equation for s, defining 



' E LPM k 
8E(E - k) ' 



(78) 



Then, £ depends only on s'\ 

=2 (s'< y/2 Sl ) 

, »v 1 , 0.08(1 - h)\l - (I - h) 2 } . f- . . 

£( s ') = 1 + h 5= f^/fT^ (v ^ Sl < 5 < X ) 

= 1 1 (*'>!)■ (79) 

where h = In s' / In (v^si). This transformation is possible because £ varies so slowly with 
s. 

Fig. 9 compares the energy weighted cross section per radiation length, X nyda/dy, 
Eq. (72), for several electron energies. As E rises, the cross section drops, with low energy 
photons suppressed the most. The number of photons with k < E(E — k)/E LPM is reduced. 

Although it reproduces the main terms of the Bethe-Heitler equation, Eq. (72) does 
not include all of the corrections that are typically used today. Without any suppression, 
Eq. (72) becomes 

^P=^fe 2 + 2[l + (l- !/ f]}^, (80) 

matching the Bethe and Heitler (1934) result. In comparison a modern bremsstrahlung 
emission cross section (Tsai, 1974) (Perl, 1994) 



g?obh 4ar 



dk 3k 



{y 2 + 2[1 + (1 - y) 2 ]} [Z 2 (F el - f) + ZF inel \ + (1 - y) { — ± Z) 



(81) 



includes several additional terms. F e i is the elastic form factor from Sec. 2, In (184/Z 1 / 3 ), 
for interactions with the atomic nucleus. F ine i is an inelastic form factor, In (1194/Z 2 / 3 ) that 
accounts for inelastic interactions with the atomic electrons. Newer suppression calculations, 
discussed later, treat inelastic scattering separately. 

The / term accounts for Coulomb corrections because the interaction takes place with 
the electron in the Coulomb field of the nucleus. For lead / = 0.33, while F e i = 3.7. The 
Coulomb correction may be incorporated into suppression calculations by adjusting the form 
factors (Sec. VII; Baier and Katkov, 1997a). 

These corrections may be accounted for with a simple assumption (Anthony et ai, 1995). 
Since q\\ <^ mc the target mass is irrelevant for suppression purposes, and the two form 
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FIG. 9. The energy weighted differential cross section for bremsstrahlung, Xonyda^PM/ 'dy for 
various electron energies in a lead target, showing how the spectral shape changes. Electrons of 
energies 10 GeV (top curve), 100 GeV, 1 TeV, 10 TeV, 100 TeV, 1 PeV and 10 PeV (bottom curve) 
are shown. The units are fractional energy per radiation length. 

factors can be lumped together. The same can be done for the Coulomb corrections. This 
is done by scaling the radiation length to include these corrections; the standard tables of 
radiation lengths (Barnett, 1996) include these factors. 

The other difference, the final 1 — y term, is problematic because of the different y 
dependence. In a semi-classical derivation (Ter-Mikaelian, 1972, pg. 18-20), this term only 
appears in the no-screening limit, small impact parameter, high momentum transfer, limit, 
where the formation length is short. So, this term should represent a part of the cross section 
that involves large momentum transfers, and so is not subject to suppression. Because it is 
only about a 2.5% correction for large Z nuclei, current experiments have limited sensitivity 
to this point. 

In the strong suppression limit, for y <C 1, the small s approximations for 0(s) and £(s) 
lead to the semi-classical scaling 



Because £(s) varies only logarithmically with s, this limit, corresponding to £(s) = 2, is only 
reached for very strong suppressions, S < 10~ 3 . As Fig. 1 shows, the normalization differs 
from the semi-classical results, but it, as previously mentioned, can depend on the how the 
absolute cross section is treated. Moreover, for a direct comparison, it might be fairer to 
use £ = 1, since the semi-classical calculations neglect the relevant non-linear scaling. This 
changes the coefficient to 3/\/2 ~ 2. For lower energies, the strong suppression limit is 




(82) 
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S ^(>-» J3M ) < 83 > 

where the recursion for £(s) has been removed, and si calculated for lead. It is clear that 
the correction for £ is significant, particularly for small S. 

These calculations contain some sources of error. Because of the Fokker-Planck method 
and the problematic mean scattering angle; these results are of only logarithmic accuracy. 
One numerical hint of inaccuracy can be seen by comparing Migdal's formula with the Bethe- 
Heitler limit. For k LPM < k < 1.3k LPM , 0(s)£(s) rises slightly above 1 and the Migdal result 
is slightly above Bethe Heitler. The maximum excess is about 3%. This can serve as a crude 
estimate of the expected accuracy of Eq. (72). 

Migdal also considered dielectric suppression. Since it only occurs for y 1, only the <fi 
term is relevant; Migdal replaced (f>(s) with (p(sT)/T, where T = 1 + k 2 /k 2 , to get 



da LPM _ lQar 2 £(s) <f>(sT) 2 ( 184 \ 



The same substitution applies for the simplified polynomials (Stanev et ai, 1982). This is 
close to the Ter-Mikaelian result, although it misses the modified form factor in Eq. (45). 

Equation (72) can also be found with a classical path integral approach (Laskin, Maz- 
manishvili and Shul'ga, 1984). The radiation for a given trajectory, Eq, (24), can be averaged 
over all possible paths, In the limit k <C E, required by the classical nature of the calcula- 
tion, this reproduces Migdal's result. Dielectric suppression can also be included in the path 
integral approach (Laskin, Mazmanishvili, Nasonov and Shul'ga, 1985). 

Migdal notes that, for thick slabs, the photon angular distribution is dominated by the 
electron multiple scattering. However, nothing in Migdal's calculation should change the 
semi-classical result that suppression should be reduced for photons with # 7 > 1/7; this 
effect may be visible in thinner targets. Pafomov (1965) discussed the energy and angular 
distribution of photons emerging from thin slabs (T <C X ), producing a complex set of 
results. 



B. Pair Creation 

The cross section for pair production may be found simply by crossing the Feynman 
diagram for bremsstrahlung, as shown in Fig. 6. Migdal used this crossing to calculate the 
cross section for pair production: 

fcLPM( y e+e " ) = ^|P(g(s-) + + (i - f ma) (85) 

where 

_ _ / E LPM k mc 2 

S ~ V 8E(k - E)m * w (86) 

and G and are as previously given. In the limit s ^> 1 there is no suppression, while s — > 
indicates large suppression, both leading to the appropriate semi-classical result. These 
equations can be simplified as was done with bremsstrahlung (Stanev, 1982). 
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FIG. 10. The differential cross section for pair production, Xonda^PM/dy in lead for various photon 
energies, showing how the spectral shape changes. Cross sections are plotted for photons of energies 
1 TeV (top curve), 10 TeV, 100 TeV, 1 PeV, 10 PeV, 100 PeV and 1 EeV (bottom curve). 



Figure 10 gives the pair production cross section in lead for a number of photon energies. 
As k increases above Elpm, the cross section drops, with asymmetric pairs increasingly 
favored. For E LPM < E{k — E)/k < 1.3E LPM , <f>(s)£(s) > 1 and the Migdal cross section 
rises above the Bethe-Heitler case, matching the bremsstrahlung overshoot. It also causes 
the 'hiccup' present at k — E LPM in Fig. 2. 



C. Surface Effects 

Migdal's calculations only apply for an infinite thickness target. Several authors have 
calculated the transition radiation due to multiple scattering, based on a Fokker-Planck 
approach. Although these calculations begin with the same approach, the details differ 
significantly, as do the final results. Gol'dman (1960) added boundaries to Migdal's Fokker- 
Planck equation, calculating the emission inside and outside the target, and an interference 
term. Outside the target, there is no emission, while inside he reproduced Migdal's result. 
For the surface terms, the photon flux for y <C 1 and s <C 1 is 

^ = 4lni^ln(^) (87) 
dk irk s 2nk k y J 

per surface. This result is simlar to the semi-classical result in Eq. (52). 

Ternovskii (1960) considered the dielectric effect as well as multiple scattering, again 
starting from Migdal's kinetic equation. He considered the entire range of target thicknesses, 
including interference between closely spaced boundaries. By comparing the radiation inside 
the target with an interference term, he found that for T < aXo/2n, multiple scattering is 
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insignificant, and the Bethe-Heitler spectrum is recovered. This is a much broader range 
than the limit T < X /1720 presented in section II. F. His results for intermediate thickness 
targets, Eq. (53), apply for aX /2n < T < aX /2ns, also a wider range than in Sec. II. F. 

For thicker targets, with T 3> aX /2irs, if the dielectric dominates (s > 1 or s(k p /k) 2 > 
1), he obtains Eq. (50). Where multiple scattering dominates (s < l,s(k p /k) 2 <C 1), then 

dN__ a fk 2 n E 2 + (E-k) 2 \ x 



dk -2.kVE 2+2 E^rfs (88) 

where '% ~ 1'. This result is similar to Eq. (87), but covers the complete range of y. 
Neglecting the logarithmic term, the spectrum matches that of unsuppressed bremsstrahlung 
for a target of thickness ~ 0.5o;Xo. This has a different E and k dependence than would 
result from 1 If of Bethe-Heitler radiation. 

Unfortunately, Eq. (88) is difficult to use. It only applies for s(k p /k) 2 <C 1; no solution 
is given for the intermediate region s(k p /k) 2 ~ 1. Also, x is poorly defined. If the equation 
is extended to s(k p /k) 2 = 1, then choosing x — 1 creates a large discontinuity. The discon- 
tinuity disappears for x — 0.42, but this also eliminates transition radiation in the region 
0.17 < s < 1. 

Garibyan (1960) extended his previous work with transition radiation to include multiple 
scattering. He found that multiple scattering dominates for k > k cr , matching the semi- 
classical result for infinite targets. If £(s) is neglected, this is the same cutoff found by 
Ternovskii. Where multiple scattering dominated, he reproduced Eq. (87). Elsewhere, he 
found the usual transition radiation from the dielectric of the medium. 

Equations (87) and (88) are negative for s > 1; common sense indicates that they should 
be cut off in this region where no transition radiation is expected. However, Pafomov (1964) 
considered this evidence that these works were wrong. He stated that they incorrectly 
separated the radiation into transition radiation and bremsstrahlung. He calculated the 
transition radiation as the difference between the emission in a solid with and without a 
gap, again starting with the Gol'dman's initial formulae. Pafomov found that multiple 
scattering dominated transition radiation for k > k cr . Surprisingly, he found that multiple 
scattering affected transition radiation even for E < ydieEhPu-, with an additional term 
dN/dk = a/nk(kLPM/k) 2 added to the conventional result when k > k p . For /c^pm > k pi 
he found different formula for different k. For k < k cr 



dN _ a _ k 
dk irk k 



n = i»7f (89) 



while for k cr < k < k^pM, 



dN a 2 k LPM , af ,, 
-dk=-k ln 3 ] J— (90) 

Counterintuitively, Pafomov predicted that for k > kiPM, there is still transition radiation, 
with 

dN_ = (^_ Sk 2 LPM 

dk irk 21k 2 ' 1 ' 
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These two equations do not match for k = kiPM- However, a numerical formula, not 
repeated here, covers the entire range smoothly. It is worth noting that Eq. (91) is quite 
close to the semi classical result for k < k LPM . 

Fig. 11 compares the transition radiation from a single surface predicted by 
Gol'dman/Garibyan, Ternovskii (using x — 1-0), and Pafomov. Except for Pafomov, these 
calculations predict radiation up to k = kiPM = E 2 /Elpm- For a flat kdN/dk spectrum, 
the total energy radiated per surface rises as E 2 , as with the semi-classical approach. For 
large enough E, electrons lose most of their energy to radiation when crossing an interface. 

However, these equations fail before this point. For high enough electron energies, these 
calculations predict that each electron should emit several photons per edge traversed. Since 
the formation lengths for the various photon emissions will overlap at the edge, this is really 
a higher order process, not yet treated properly by calculations. 

Unfortunately, these calculations appear to do a poor job of fitting the data. Fig. 19 
shows that they predict transition radiation considerably above the data from SLAC E-146. 



IV. BLANKENBECLER AND DRELL FORMULATION 

Blankenbecler and Drell (1996) calculated the magnitude of LPM suppression with an 
eikonal formalism used to study scattering from an extended target. The approach was 
originally developed to study beamstrahlung. One major advantage of their approach is that 
it naturally accommodates finite thickness slabs, automatically including surface terms. 

They begin with a wave packet that scatters while moving through a random medium. 
For each electron path, they calculate the radiation, based on the acceleration of the elec- 
tron. The radiation is calculated for all possible paths and averaged. One clear conceptual 
advantage of this calculation is that it does not single out a single hard scatter as causing 
the bremsstrahlung; instead all of the scatters are created equal. This differs from Lan- 
dau and Pomeranchuk and Migdal, who found the rate of hard scatters which produced 
bremsstrahlung, and then calculated the multiple scattering in the region around the hard 
scatter to determine the suppression. For a thick target, in the strong suppression limit 
s<l, they predict that the radiation is ^3^/8 (about 8%) lar ger that of Migdal. 

One advantage of this calculation is that it treats finite thickness media properly. There 
are three relevant length scales: lf (Eq. (5)), T and aX . Blankenbecler and Drell (1996) 
called aX the mean free path for elastic scattering, based on counting vertices in Feynman 
diagrams. This correspondence does not stand up to more detailed examination. However, 
the arguments in the paper are not affected. These variables are combined into two ratios: 

and 

Tbd = (93) 

where T B d is the target thickness, in mean free paths, while N B d is the number of formation 
lengths per mean free path; when N B d is large suppression is strong, while N B d < 1 
corresponds to a single interaction per electron, i.e. the Bethe Heitler regime. 
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FIG. 11. Comparison of the edge effects predicted by Ternovskii (with x = 1) (solid line) and 
Pafomov (long dashes), for a single independent edge, compared with conventional transition ra- 
diation (short dashes), (a) shows the radiation from a 100 GeV electron in lead. The jump in the 
Ternovskii curve at k = 350 keV is at s{k p /k) 2 ) = 1, while the drop in the Pafomov curve at 630 
keV is at k = k cr . Panel (b) shows the Pafomov and Ternovskii predictions for a 10 TeV electron 
in lead; both curves are much smoother. Conventional transition radiation is negligible in (b). 
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The eikonal approach finds the wave function phase and momentum difference between 
different points on the electron path. These differences are then used to find the radiation 
for that length scale. The probability of emission from the target is 

Here, k± = /ccos(# 7 ) is the photon perpendicular energy, S(z2,zi) is a sum over electron 
polarization (spin flip and no spin flip) of the squared perpendicular momentum acquired 
by the electron due to scattering, and d^/dz is the differential phase difference due to the 
scattering. The d?k± integration is equivalent to an integration over photon angles. In the 
absence of scattering, d§/dz is constant, and then the integrals must be carefully evalu- 
ated; Blankenbecler and Drell introduce a convergence factor to insure that the boundary 
conditions are correct. 

Because the integrals include all possible values of z, three regions must be considered: 
before the target (denoted +), inside the target (denoted 0), and after the target (— ). For 

the double integral, there are 9 possible combinations, of which two (++ and ) are clearly 

zero. Time ordering eliminates ( — h), (—0) and (0+), leaving the bulk term (00), two single 
surface terms (+0) and (0—) and an interference term (H — ). For thick targets with T BD 3> 1, 
the interference term vanishes, and the bulk term dominates over the single surface terms. In 
this case, their calculations reproduce the Migdal results, with the 8% higher cross section. 
For thinner targets, the surface terms are more important. Figure 12 compares these terms, 
using a suppression form factor S = F(N BD ,T BD ) relative to Bethe-Heitler. 

For thin targets, T B d < 1, where N BD > 1, the interference term (H — ) dominates, 
demonstrating a large transition radiation. N BB ./T BB < 1 reduces to the Bethe Heitler free 
particle case. 

For a thick target, T BB) > 1, the central interaction (00) region is important. For large 
N BD) their results are similar to those of Migdal. For very large N BD) with N BD > T BD) the 
formation length is longer than the target, and both the (00) and mixed regions contribute. 
Fig. 13 compares the suppression as a function of N BD (which is proportional to 1/k), for 
a variety of T. For large T BB , suppression increases with N BB . When T BB decreases, 
suppression drops, eventually becoming almost independent of N BB , at a value similar to 
that predicted by Shul'ga and Fomin (1996). 

There are a few caveats in this calculation. The eikonal approach assumes that the 
potential is smooth enough. Blankenbecler and Drell used a Gaussian scattering potential, 
which underestimates the rate of large angle scattering. Second, they assumed that the 
wave function phase and amplitude fluctuate independently. Blankenbecler (1997b) showed 
that the correlation between amplitude and phase reduces the emission by a further 5-15%. 
Figure 13 compares curves with and without the correlation. The calculation has also been 
extended to include multiple slabs separated by a gap (Blankenbecler, 1997a). 

Calculating the emission from a slab as a whole is problematic, These results assume that 
there is either zero or one interaction per incident electron. But, for a typical bremsstrahlung 
cross section of 10 photons per X , the relative probability of getting 2 interactions compared 
with 1 interaction is ~ 20T/X , so single interactions are only prevalent for slabs with 
T < 0.05X . This problem makes it difficult to apply these results in many real-world 
situations. 
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FIG. 12. The Blankenbecler and Drell form factor F(N B d,T B d) for a target with T BD = 10, 
showing the contributions for the bulk emission (00) and the transition and interference terms 
(0— ), (+0) and (H — ). For a fixed target thickness, Nbd ~ V/o ~ k/E 2 ; the edge effects are 
largest for small Nbd corresponding to small k. From Blankenbecler and Drell (1996). 
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FIG. 13. The Blankenbecler and Drell form factor F(N BD ,T BD ) for 25 GeV electrons in 3 thick- 
nesses of gold targets. F(lpm), assumes that the wave function amplitude and phase fluctuate 
independently, while F(tot) includes the correlation between amplitude and phase, further reduc- 
ing the emission. From Blankenbecler (1997b). 
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V. ZAKHAROV CALCULATION 



Zakharov (1996a,b, 1998a,b) transformed the problem of multiple scattering into a two 
dimensional (impact parameter and depth in the target) Schroedingers equation. The imag- 
inary potential is proportional to the cross section for an e + e~ pair scattering off of the 
atom, which is itself simply related to the bremsstrahlung cross section. This equation was 
solved using a transverse Greens function based on a path integral. This impact parameter 
approach is complementary to Migdal's momentum space approach. Because this approach 
allows for arbitrary density profiles, it naturally accommodates finite thickness targets. 

Zakharov (1996a) calculated the radiation due to a simple Coulomb potential. The calcu- 
lation is keyed to the scattering cross section for a e + e~ dipole of separation p, a(p) = C(p)p 2 . 
In the strong suppression limit, C varies slowly with p. Zakharov then found the frequency 
Q for a harmonic oscillator in a potential which would reproduce this scattering cross sec- 
tion. This is roughly equivalent to describing the multiple scattering with a Gaussian. The 
radiation is governed by a parameter 77 = Qlf. For an infinitely thick target, the results are 
almost identical to Migdal (Zakharov, 1998a); functions of 77 that match Migdal's <f>(s) and 
G(s), for 77 = 1/ \^8s. The only difference is the slowly varying part of the cross section: £(s) 
for Migdal and C \y p e s) / C '{Tic / m) for Zakharov; p e s is the impact parameter where radiation 
is largest. In the Bethe- Heftier limit, p c s ~ Hc/ym. 

In the limit y — > 0, for strong suppression 



where a is the Thomas- Fermi screening radius Numerically, p e ff( m ) ~ 10 9 /[E(eV)y 3 ] 1 ^ for 
lead. These equations are valid for 200/Z 1 / 3 < log(2a/yp eS ) < 1-5 x 10 5 /Z 2 / 3 ; corresponding 
to Migdal's f (s) = 1 and f (s) = 2. 

Except for the '1' in £(s), this equation has the same form as Migdal's Eq. (83). Although 
unimportant to the theory, the '1' greatly reduces the effect of the slowly varying term. 
Numerically, 2a/yp cS = 6.8 x 1CT 7 E(eV)/y for lead; other solids aren't too different. For a 
1 TeV electron beam, as y varies from 10~ 4 (a typical y<jj e ) to 1CT 2 (an arbitrary upper limit 
to 'low y'), the Zakharov logarithm varies by about 20%, while the change in £(s) is much 
smaller. This variation should be measureable. 

For moderate suppression, a more detailed treatment is required. Because C(p) varies 
more quickly with p, the harmonic oscillator approximation fails and the actual potential 
must be used (Zakharov, 1996b). Zakharov used separate screened elastic (a = 0.83O0.Z -1 / 3 ) 
and inelastic (a = 5.2a Z~ 2 / 3 ) potentials, reproducing the appropriate unsuppressed 
bremsstrahlung cross sections. The separate potentials are most important for low Z nuclei. 
However, because of the small recoil, the separate form factors have a small effect on sup- 
pression. The more complicated potential requires additional integrations. Because of this, 
this approach can so far be used only for finite thickness targets. 



da z _ 2a 2 Z 2hc\og{2a/yp eS ) 



(95) 




where 




(96) 
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FIG. 14. Suppression factor S(h,k) found by Zakharov. The curves are for a 25 GeV electron, 
with h = T/lfQ. The h = 10 results are very close to the infinite thickness limit. The dip around 
h = 1 for the highest photon energies is due to interference between the two target surfaces. From 
Zakharov (1996b). 

Figure 14 shows Zakharov's suppression factors for finite target thicknesses in a 25 GeV 
electron beam. This demonstrates how suppression increases with target thickness. The 
thickest target, with h = 10 is very close to the infinite thickness limit. Zakharov (1998b) 
added a correction to allow for multiple photon emission, and did a detailed comparison 
with SLAC E-146 data. Except for the carbon targets in 25 GeV electron beams, he found 
good agreement with the data for photon energies k > 5 MeV (above the region of dielectric 
suppression) . 

Zakharov (1997b) presents a few results for multiple-slab configurations, finding a smaller 
interference term than Blankenbecler (1997a). The two results would agree better if Blanken- 
becler included the amplitude-phase correlation in his multiple slab calculations. 

Zakharov's results for gluon bremsstrahlung from a quark will be discussed in Sec. XII. 



VI. BDMS CALCULATION 

The BDMS group (R. Baier et ai, 1996) started with the Coulomb field of a large number 
of scatterers, with the atomic screening modelled with a Coulomb potential cut off with a 
Debye screening mass fi. An eikonal approach is used to account for the large number of 
scatters. 

The critical variable in this calculation is the (dimensionless) phase difference between 
neighboring centers, 

X fi 2 c 3 k 

= ^m&r (97) 
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For QED calculations, \i = TiZ 1 ^ /a^c reproduces Thomas- Fermi screening. The authors also 
define a coherence number v = alfo/X , the number of scatters required for the accumulated 
phase shift to equal 1. This v is similar to Blankenbecler and Drells Nbd- A large phase shift 
between interactions, k > 1, corresponds to the Bethe-Heitler limit. Because this approach 
assumes massless electrons, its 'Bethe-Heitler limit' does not completely match the usual 
Bethe-Heitler formula. 

In the factorization limit, T < lf , where the entire target reacts coherently, their results 
match Ternovskii (1960): 

where q tot is the total perpendicular momentum acquired by the particle while traversing 
the target due to multiple scattering. Here, the mass m is introduced to remove a collinear 
divergence. 

In the LPM regime, (k <C 1), the bremsstrahlung cross section is similar to Eq. (95). The 
logarithmic term is due to non-Gaussian large angle (9 > 1/7) Coulomb scatters. Because 
of these large scatters, the mean squared momentum transfer is poorly defined, and, in fact, 
diverges logarithmically; this logarithm appears in the suppression formula. 

da a 



dk nkX 



Kln(l//c). (99) 



They comment that they neglect a logarithmic factor under the logarithm; if the corre- 
sponding term is removed from Zakharov's formula, the two results have the same functional 
dependence. If one identifies Elpm = Kfi 2 /2, as their paper indicates, then the radiation 
takes the form 

do a I kE LPM ~W 

^~^V~^ l0g *w (100) 



VII. BAIER 81 KATKOV 

Baier and Katkov (1997a) also studied suppression due to multiple scattering, trying to 
reach an accuracy of a few percent, by including several corrections omitted previously. They 
begin with the scattering from a screened Coulomb potential, in the same impact parameter 
space used by Zakharov. Coulomb corrections are included to account for the motion of 
the screening electrons, along with separate potentials for elastic and inelastic scattering. 
Finally, they allow for a nuclear form factor, with an appropriately modified potential for 
impact parameters smaller than the nuclear radius. 

They find the electron propagator for a screened Coulomb potential, in the Born ap- 
proximation. This assumes Gaussian distributed scattering, and reproduces Migdal's result. 
The electron propagator is then expanded perturbatively, with a correction term which ac- 
counts for both large angle scatters and Coulomb corrections to the potential. Without the 
Coulomb corrections, the first order result is similar to that of Zakharov (V. Baier, 1998). 

Coulomb corrections are incorporated by adjusting the parameters in the potential. With 
Coulomb corrections, the screening radius becomes a 2 = 0.81Z~ 1//3 a exp[0.5 — f(Za)]. while 
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the characteristic scattering angle changes from 9\ = he/ Ea^ to 62 = hc/Eao exp[f(Za) — 
0.5)]. Here, f(Za) is the standard Coulomb correction (Tsai, 1974). This scaling accounts 
for the extra terms in Eq. (81) compared with Eq. (80). For heavy nuclei, a 2 is about 20% 
larger than the Thomas-Fermi screening radius. 

Higher order terms may also be calculated for the electron propagator. The ratio of the 
first two terms of the expansion is 0.451 divided by a logarithmic term, so the series should 
converge reasonably rapidly. 

Inelastic scattering can be included by adding a term to the scattering potential, changing 
the charge coupling from Z 2 to Z 2 + Z and further modifying the characteristic angle, to 
9 e = 9 1 exp [Z/(l + Z){f{aZ) - 1.88) - 0.5]. 

Dielectric suppression is included by modifying the potential, with a replacement similar 
to Migdal's. For k <C k p , the spectrum is similar to Ter-Mikaelian's, but with power law 
and Coulomb corrections: 

ln (i^) + n-'H < 101 > 

where f(Zct) is the standard Coulomb correction (Tsai, 1974). 

Baier and Katkov considered the extremely strong suppression regime, neglected by 
Migdal, where the finite nuclear radius becomes important. This region is reached at the 
lowest E for ?/ < 1, where dielectric suppression otherwise dominates. There, limiting q 
(i.e. q\\) to H/Ra changes the form factor for k < kpR^Xe- In lead, this corresponds to 
y < 2 x 10~ 6 . Then, 

For y = 10~ 6 in lead, this is about 25% larger than Eq. (101). This is probably measurable, 
although transition radiation and backgrounds will be very large. 

Baier and Katkov then considered targets with finite thicknesses, breaking down the 
possibilities in a manner similar to Blankenbecler and Drell, with a similar double integral. 
For relatively thick targets, T > lf , the results are consistent with Ternovskii (1960), but 
with additional terms for the Coulomb correction. For T < lf and strong suppression, 

dN a\k 2 ( (E — k) 2 \ ( 1 , n . t . mn . 1 0.578\1 , . 

lk = -kW + V + L -r L ) ( (1 + 2A^**> ~ °' 5781 + TA ~ 1 + — )] (103) 

where A = nZ 2 a 2 nTh 2 /m 2 c 2 (L t + l-2 x 0.578) and L t ~ 2 In (2a 2 /\ e pt), with p t the (scaled) 
minimum impact parameter that contributes to the form factor integral. Unfortunately, this 
equation diverges as T — > 0. 

V. Baier and Katkov (1997a) compared their photon spectrum with SLAC E-146 data 
from a 0.02X thick tungsten target in 25 and 8 GeV electron beams, and find good agree- 
ment. However, a target this thick has a substantial (roughly 20%) correction to the photon 
spectrum to account for electrons that undergo two independent bremsstrahlung emissions. 
This correction is not included in their calculation, and the good agreement with data is sur- 
prising. A later work (Baier and Katkov, 1999) included multiphoton effects, and also finds 
good agreement with the data; the difference between the two calculations is not explained. 
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da 
dk 



IbZ^ar^k 
3^2 



(102) 



V. Baier and Katkov (1997b) considered thinner targets with T ~ IfQ. They compare 
their calculations with SLAC E-146 data, for 0.7 % X Q thick gold targets in 8 and 25 GeV 
beams, and also find good agreement. Because this target is much thinner, the multiple 
interaction probability is greatly reduced and the agreement is not surprising. 

VIII. THEORETICAL CONCLUSIONS 
A. Comparison of Different Calculations 

The calculations discussed here used a variety of approaches to solve a very difficult 
problem. Because the underlying techniques are so different, it is difficult to compare the 
calculations themselves. However, some general remarks are in order. 

All of the post-Migdal calculations are done for a finite thickness slab, integrating both 
bulk emission and transition radiation. Unfortunately, the finite thickness slab calculations 
are not easily usable by experimenters, because they assume each electron undergoes at 
most one interaction in the target. Multiple interactions are easily accounted for in a Monte 
Carlo simulation. However, the simulation must be able to localize the photon emission, 
while these calculations are for the slab as a whole. Because of the edge terms, it isn't 
correct to simply spread the emission evenly through the slab. This limits their direct 
applicability to thin slabs. 

The newest calculations by Zakharov and Baier and Katkov include multiple photon 
emission. However, these calculations are still for bulk targets, and are not very amenable 
to complex geometries. For example, these calculations could not be used to model an 
electromagnetic shower. 

Here, we will use the Migdal approach as a standard for comparison. Although Za- 
kharov's approach is very different from Migdal, he reproduces Migdal's result in the appro- 
priate limit. However, Zakharov incorporates some further refinements which should lead 
to increased accuracy. Baier and Katkov have a similar approach to Zakharov, and include 
a number of additional refinements, especially for bremsstrahlung at very low y. 

The other works have very different genealogies. Because the BDMS result does not 
reproduce the Bethe- Heftier limit, it is more relevant to QCD than QED. However, if the 
formalism is extended to include hard radiation, and some additional graphs are added, then 
BDMS becomes equivalent to the Zakharov formalism (R. Baier et al, 1998b). 

For Blankenbecler and Drell, the only obvious point of comparison is the potential; 
Zakharov (1996b) states that the Blankenbecler and Drell potential does not match the 
Coulomb potential, and will not show the logarithmic dependence given by the slow variation 
ofC(p) or£(s). 

The results can also be compared numerically. Unfortunately, these calculations are 
very complex and the descriptions frequently lack adequate information for independent 
computation. So, it is necessary to rely on the results given by the authors. One point 
of comparison is defined by SLAC E-146 data on a thin target: 25 GeV electrons passing 
through a 0.7 % Xo thick gold target. The E-146 data agreed well with Migdal as long as T > 
If] at lower k, Ternovskii's Eq. (53) matched the data. Blankenbecler and Drell, Zakharov, 
and Baier and Katkov all showed good agreement with this data. Zakharov initially added 
in a 7% normalization factor to match the data. Some correction is required, because, 



42 



even for a 0.7% Xq target, the multiphoton pileup 'correction' is still several percent. With 
a correction for multiphoton emission Zakharov (1998b) found normalization coefficients 
average 1, except for the uranium targets. Zakharov also changed his approach to Coulomb 
corrections between the two works. The other authors do not discuss normalization. Overall, 
in this energy range, the different approaches agree with each other to within about 5%. 

Because of the different logarithmic treatment, Migdal, Blankenbecler and Drell, and 
Zakharov will scale slightly differently with energy. Changes depending on target thickness 
are more complicated; the surprising agreement found by Baier and Katkov for the 2% X 
tungsten data, where multiple interactions are a 20% correction, shows that there there 
are significant uncertainties in scaling the results with target thickness. It would also be 
interesting to compare the calculations for a low Z target, where the E-146 data showed 
some disagreement with Migdal. 

B. Very Large Suppression 

One weakness of all of these calculations is that they only consider the lowest order dia- 
grams. For fixed y, If rises with E, even with LPM suppression. At high enough energies, 
the formation zones from different emissions will overlap, and any lowest order calculation 
will fail. Dielectric suppression is strong enough that // decreases and localization improves 
with increasing suppression, so it is less subject to this problem. However, for multiple 
scattering, a method of dealing with higher order terms is needed. While the radiative 
corrections to bremsstrahlung are known (Fomin, 1958), they were not computed with sup- 
pression mechanisms in mind. 

Even neglecting the overlapping formation zones, when suppression is large, higher order 
terms are important, because higher order processes involve larger q\\, and hence are less 
subject to suppression; Sec. II. H illustrated this for direct pair production. So, when 
suppression is strong (S ~ a), current calculations are suspect. For QCD calculations, 
discussed in Sec. XII the problem is much worse because of the large coupling constant, a s . 

It is worth noting that suppression can affect other processes, via radiative corrections. 
The total cross section for Coulomb scattering, for example, is the sum of the elastic cross 
section plus the cross section for bremsstrahlung where the photon is not observable. Since 
suppression can affect the latter process, it can indirectly change the elastic scattering cross 
section (Ter-Mikaelian, 1972, pg. 135). Inelastic processes, which leave the target atom in 
an excited state can also contribute significantly to the cross section when suppression of 
the lowest order diagrams is large. 

C. Classical and Quantum Mechanical Approaches 

There has been some controversy as to how well classical electrodynamics can predict 
LPM suppression. The original Landau and Pomeranchuk (1953a,b) calculations were com- 
pletely classical; they failed at the same point that classical bremsstrahlung calculation fail: 
when k ~ E. The path integral approach (Laskin, Mazmanishvili and Shul'ga, 1984) goes 
further, and reproduces Migdal's result in the same limit. Of course, at least semi-classical 
calculations are necessary for pair production. 
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The main advantage of quantum mechanical calculations is that it covers the complete 
range of y. The newest calculations include features not found in the classical approaches. 
However, some of the advances, such as separate elastic and inelastic potentials could cer- 
tainly be included in a classical calculation. 

In the future, quantum approaches appear necessary to calculate the higher order terms 
that become important at extremely high energies, when formation lengths from different 
emissions begin overlapping. Unfortunately, these calculations seem extremely difficult, and 
are probably a long way off. 

However, the classical approach to LPM suppression was not universally accepted. Bell 
(1958) stated that Migdal's predictions conflicted with classical electrodynamics, and con- 
cluded that "any real effect of this kind is of essentially quantal origin." He pointed out 
that the classical radiation, Eq. (6), is positive definite, and hence monotonically increasing 
for increasing pathlength. However, Bell neglected to account for the fact that, in a dense 
solid, the electron trajectory r(t) changes, and one cannot simply sum the radiation due 
to interactions with separate nuclei. Even classically, emission from the different pieces of 
electron pathlength can interfere, in a manner similar to other classical interference effects. 

IX. EXPERIMENTAL RESULTS 

Bremsstrahlung or pair creation suppression can be studied with high energy electron or 
photon beams. Because pair creation suppression requires photons with k > E LPM , beyond 
the reach of current accelerators, pair creation has been studied only with cosmic rays, with 
consequently very limited statistics. The best suppression studies have used electron beams 
at accelerators. Besides the LPM effect, these beams have been used to study dielectric 
suppression. 

A. Cosmic Ray Experiments 

The first tests of LPM suppression came shortly after Migdal's paper appeared. These 
experiments used high energy (k > 1 TeV) photons in cosmic rays, and studied the depth 
of pair conversion in a dense target. The earliest experiments looked for a deficiency of low 
energy electrons in electromagnetic cascades (Miesowicz, Stanisz and Wolther, 1957). This 
analysis is difficult because of uncertainties in the incoming photon spectrum and electron 
energy measurement. 

Fowler, Perkins and Pinkau (1959) studied isolated pair production in an emulsion stack. 
They measured the shower energy and the distance between the initial conversion and the 
first daughter pair resulting from the primary. They observed 47 showers with k > 1 TeV. As 
k increased, the distance between the conversions rose, as predicted by Migdal, contrasting 
with with decrease predicted by Bethe- Heitler. 

Varfolomeev and collaborators (1960) used a similar approach, albeit on a smaller data 
set, with similar results. They selected events where the energy of the daughter pair was 
very small (y <C 1), and found qualitative evidence for suppression; they did not differentiate 
between LPM and dielectric suppression. 

Lohrmann (1961), studied emulsion exposed on high altitude balloon flights. Compared 
to earlier results, Lohrmann studied lower energies, where the LPM and Bethe-Heitler cross 
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sections are much closer. However, the increased photon flux allowed for much better statis- 
tics, making up for the smaller cross section difference. Lohrmann used a variety of analysis 
techniques, all of which supported Migdal's results. More recently, long duration balloon 
experiments have gathered somewhat larger data samples, up to 120 events, with similar 
results (Strausz et al, 1991). 

Kasahara (1985) studied the development of ~100 TeV showers in lead/emulsion cham- 
bers exposed to cosmic rays on Mt. Fuji. His analysis was based on the fact that elec- 
tromagnetic LPM showers with E ^> E LPM penetrate considerably more material than 
Bethe-Heitler showers. Kasahara found that the shower development profiles matched simu- 
lations based on Migdal, but differed from his Bethe-Heitler simulations. The study suffered 
from a background from hadronic and/or multiple photon showers, both also subject to 
large fluctuations. Kasahara removed 5 obvious background events from his initial 19 event 
sample, but warned that some of the remainder could also be hadronic. This study is of 
special interest as it is the only study yet in the region E 3> E LPM , where suppression affects 
shower development. 

All of the air shower experiments suffered from some common problems, by far the largest 
being the poor statistics; high enough energy photons are not common. In addition, because 
of the limited emulsion thickness, they did not consider surface effects. Finally, uncertainties 
in the photon spectrum complicated the analysis. For these reasons, these experiments are 
at best qualitative verifications of LPM suppression. 

B. Early Accelerator Experiments 

By the 1970's, 40 GeV electron beams were available at Serpukhov, and the first 
accelerator based test of LPM suppression was done there (Varfolomeev et al, 1975). 
Bremsstrahlung photons from 40 GeV electron beams striking dense targets were detected 
in an sodium iodide calorimeter. The electrons were magnetically bent away from the 
calorimeter. Photons with 20 MeV < k < 70 MeV emitted from carbon, aluminum, lead 
and tungsten targets were studied. Below 20 MeV, synchrotron radiation from the bending 
magnets dominated the measurement, while above 70 MeV, the experiment was insensitive 
to the spectral change predicted by Migdal. 

The experiment suffered from several limitations. Because the electron beamline was 
in air, and included several scintillation counters used as triggers, there was a significant 
light element bremsstrahlung background. The experimenters also mention a significant 
background due to muon contamination in their beam. 

The collaboration presented their data in terms of ratios of photon spectra: 
lead/aluminum and tungsten/carbon. This may have been done to account for events con- 
taining several bremsstrahlung photons from a single electron. The data showed suppression 
in the region that Migdal predicted. However, the degree of suppression was larger than 
Migdal predicted, although within the large errors. 

CERN NA-43 was an experiment dedicated to studying channeling radiation from elec- 
trons and positrons in crystals (Bak et al, 1988). In channeling, electrons or positrons 
travel along the crystal rows, and hence are strongly affected by the coherently adding 
fields from the atom rows. At large angles to the axes, the coherence disappears, and nor- 
mal bremsstrahlung occurs. In this large angle regime, they saw suppression which they 
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attributed to LPM suppression. They also observed suppression consistent with multiple 
scattering for electrons incident along one of the crystal rows. This is slightly different from 
suppression in an amorphous material; here the suppression is really the loss of coherent 
enhancement. 

The first experimental studies of dielectric suppression were part of a larger study of 
transition radiation (Arutyunyan, Nazaryan and Frangyan, 1971). Although the experiment 
focused on studies of emission from stacks of thin radiators, the experimenters also measured 
the radiation from 0.25 and 2.8 GeV electrons traversing relatively thick glass and aluminum 
targets. Although few details were given, they appeared to observe suppression in the 
expected energy range. 

C. SLAC E-146 

In 1992, the E-146 collaboration at Stanford Linear Accelerator Center (SLAC) proposed 
an experiment to perform a precision measurement of LPM suppression, and to study di- 
electric suppression (Cavalli-Sforza et al, 1992). The experiment was conceptually similar 
to the Serpukhov experiment, but heavily optimized to minimize background. To minimize 
the statistical errors, the experimenters collected a very large data set. The experiment was 
approved in December, 1992, and took data in March-April, 1993. 

1. Experimental Setup 

Figure ?? shows a diagram of the experiment. An 8 or 25 GeV electron beam entered 
SLAC End Station A and passed through a thin target. The targets used are listed in Table 
III. The beam was then bent downward by a 3.25 T-m dipole magnet, through six wire 
chamber planes which measured electron momenta and into an array of lead glass blocks 
which accurately counted electrons. Produced photons continued downstream 50 meters 
into a BGO calorimeter array. To minimize backgrounds, the electron path visible to the 
calorimeter and the photon flight path were kept in vacuum. 
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FIG. 15. Diagram of the SLAC-E-146 apparatus. From Anthony et al. (1995). 



TABLE III. The targets used in SLAC E-146. These &lpm values are in MeV, and are half those 
used by the E-146 collaboration. 



Target Material 


Thickness (Xo) 


k LPM (25 GeV) 


k LPM (8 GeV) 


Carbon 


2%, 6% 


4.3 


0.4 


Aluminum 


3%, 6% 


7.8 


0.8 


Iron 


3%, 6% 


48 


4.8 


Tungsten 


2%, 6% 


236 


24 


Gold 


0.07%, 0.7%, 6% 


250 


26 


Lead 


2% 


148 


15 


Uranium 


3%, 5% 


236 


24 



Blank 



The calorimeter comprised 45 BGO crystals in a 7 by 7 array with missing corners; each 
crystal was 2 cm square by 20 cm (18 X ) deep. This segmentation provided excellent spatial 
resolution for separating synchrotron radiation from bremsstrahlung photons. Scintillation 
light from each crystal was measured separately by a photomultipler tube (PMT). The light 
yield was about 1 detected photoelectron per 30 keV, providing good statistics down to 200 
keV. The calorimeter resolution was about 8% (FWHM) at 100 MeV, with a nonlinearity 
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of about 3%. The calorimeter temperature was monitored throughout the experiment, and 
the data corrected using the measured drifts. 

The collaboration used several methods to calibrate the calorimeter, to obtain both an 
absolute energy calibration and a crystal to crystal intercalibration. The primary tools for 
measuring the relative gain were nearly vertical cosmic-ray muons, selected by a plastic 
scintillator paddle trigger; the gain in each crystal channel was adjusted to produce equal 
signals. 

The absolute energy scale of the calorimeter was primarily determined using a 500 MeV 
electron beam. This calibration was checked with data by comparing the electron energy 
loss, measured by the wire chamber, with the calorimeter energy measurement. Because of 
the steeply falling photon spectrum and the non-Gaussian errors in the momentum mea- 
surement, this was useful only as a cross-check. 

Because of the large bremsstrahlung cross section, the experiment required a beam in- 
tensity of about 1 electron per pulse. Because it would have been very uneconomical to use 
the SLAC linac to produce a single electron per pulse, the collaboration developed a method 
to run parasitically during SLAC linear collider (SLC) operations, by using the off-axis elec- 
trons and positrons that are removed from the beam by scrapers in the beam switchyard 
(Cavalli-Sforza et al, 1994). 

Normally, about 10% of the SLC beam is scraped away by collimators in the last 200 
meters of the linac. The collimators are only 2.2 X thick, so a usable flux of high en- 
ergy photons emerged from their back and sides. Some of these photons travel down the 
beampipe, past the magnets that bend the electrons and positrons into the SLC arcs, and 
into the beam switchyard, where a 0.7 Xo target converted them into e + e" pairs. Some of 
the produced electrons were captured by the transport line optics, collimated, selected for 
energy, and transported into the end station. 

The beam worked well, with the size, divergence and yield matching simulations. At 
8 and 25 GeV, the beam intensity was adjustable up to about 100 electrons/pulse. At 1 
electron/pulse, the beam emittance was limited by the optics, with a typical momentum 
bite of Ap/p < 0.5%. The beam optics were adjusted to minimize the photon spot size at 
the calorimeter; spot sizes there were typically a few mm in diameter. The beam spot was 
stable enough and small enough that beam motion was not a major source of error. 

Data was collected and written to tape on every beam pulse (120 Hz). With the beam 
intensity averaging 1 electron/pulse, over 500,000 single electron events could be collected 
per eight hour shift. The experiment ran for a month, and good statistics were obtained 
with a variety of targets. 

2. Data Analysis and Results 

The E-146 analysis selected events containing a single electron, as counted by the lead 
glass blocks. The photon energy was found by summing the energies of hit BGO crystals 
using a cluster-finding algorithm. The cluster finding reduced the noise level by eliminating 
random noise hits. 

The experiment studied the photon energy range from 200 keV to 500 MeV, a 2500:1 
dynamic range. This exceeded the linear dynamic range of the PMTs and electronics, so 
data was taken with two different PMT gains, varied by changing the PMT high voltage. 
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The high-A; running corresponded to 1 ADC count per 100 keV and the low- A; running was 
1 ADC count per 13 keV, with the relative gains calibrated with the cosmic ray data. The 
high-A; data was used for 5 MeV < k < 500 MeV while the low- A; data covered 200 keV 
< k < 20 MeV, with a weighted average used in the overlap region. 

These two sets of data differed in several significant ways. There were large differences in 
calorimeter behavior and background levels, as well as the physics topics. For k > 10 MeV, 
photons largely interacted by pair conversion, producing an electromagnetic shower. Showers 
typically deposited energy in 3-20 crystals in the calorimeter. For k < 2 MeV, the photons 
interacted by single or multiple Compton scattering. Usually, Compton scattering deposited 
energy in a single calorimeter crystal. Sometimes, the photon Compton scattered once and 
then escaped from the calorimeter, taking some energy with it. This added a low energy 
tail to the energy deposition curve. While the high-A; data had very low backgrounds, the 
low-A; data had significant backgrounds due to synchrotron radiation, at least for the 25 
GeV beams. Finally, the two data sets emphasize different physics, with the high-A; data 
most relevant for LPM suppression, with the low-A; data more useful for studying dielectric 
suppression. For these reasons, the two sets of data were analyzed quite independently, and 
combined in the final histograms. 

3. Backgrounds and Monte Carlo 

One advance introduced by E-146 was the use of a detailed, high statistics Monte Carlo. 
The main purpose of the Monte Carlo was to understand multi-photon pileup. This occurred 
when a single electron passing through the target interacted twice, radiating two photons. 
The Monte Carlo also simulated photon absorption in the target (Anthony et al, 1997) and 
modelled the detector. Transition radiation was treated as an integral part of the physics, 
rather than a background. In addition to conventional transition radiation, the predictions 
of Ternovskii and Pafomov were included as options. 

The Monte Carlo tracked electrons through the material in small steps, allowing for the 
possibility of bremsstrahlung in the material and transition radiation at each edge. LPM 
suppression was implemented using simple formulae (Stanev et al, 1982) with dielectric 
suppression incorporated using Eq. (84). For consistency, the Bethe- Heftier cross sections 
were included by turning LPM suppression off from Migdal's formulae, rather than using a 
more modern formula. 

The expected and measured backgrounds were both small. The major background was 
synchrotron radiation from the spectrometer magnet. Synchrotron radiation was significant 
for k < 1 MeV in the 25 GeV data. Because the magnet bent the beam downward, the 
synchrotron radiation painted a stripe in the calorimeter, extending downward from the 
center. Because of the large lever arm, and because the bending started in the fringe field 
of the magnet, where the field was low, the synchrotron radiation was small near the center 
of the calorimeter. Because of the good spatial resolution of the calorimeter, most of the 
synchrotron radiation was removed with a cut on the photon position in the calorimeter; 
photons within 45° of a line downward from the calorimeter midpoint were removed; this 
cut removed most of the synchrotron radiation, along with 25% of the signal. 

Non-target related backgrounds were measured with target empty runs. The backgrounds 
were typically 0.001 photons with k > 200 keV per electron, much less than the ~ 10T/X 
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bremsstrahlung photons with k > 200 keV per electron. Target related backgrounds were 
expected to be small; photonuclear interaction rates are small, and the events are unlikely 
to appear in the E-146 analysis. 

4. Results 

Because of the high statistics and low background, the E-146 data allowed for detailed 
tests of the theory; photon spectra could be easily compared with different predictions. 

Figures 16-20 show a sampling of E-146 results. Photon energies were histogrammed 
logarithmically, using 25 bins per decade of energy, so that each bin had a fractional width 
A A; ~ 10%k. The logarithmic scale is needed to cover the 2500:1 energy range. The 
logarithmic binning dN/d\nk = kdN/dk was chosen so that a 1/k Bethe-Heitler (BH) 
spectrum will have an equal number of events in each bin, simplifying the presentation and 
statistical analysis. 

Figure 16 shows the photon spectrum (points with error bars) from 8 and 25 GeV elec- 
trons passing through 2% and 6% X carbon targets. Also shown are 3 Monte Carlo (MC) 
histograms. The top histogram (dashed line) is a simulation of BH bremsstrahlung plus 
conventional transition radiation; the transition radiation is substantial below k p , 1.4 (0.4) 
MeV at 25 (8) GeV. Above k p , the spectrum is sloped because there is a finite probability of 
a single electron interacting twice while passing through the target. Because the calorimeter 
cannot separate single photons from multiple hits, but instead measures total energy depo- 
sition, this depletes the low energy end of the spectrum (shown here), while increasing the 
number of calorimeter overflows. In the absence of multiple interactions, the Bethe-Heitler 
spectrum would be flat at (l/X )dN/d\nk = 4/3 In (k max /k m i n ) = 0.129 for bins with log- 
arithmic widths k max /k m i n = 10 1//25 = 1.096. The bin heights directly scale with the bin 
fractional width Ak/k. 

The dotted histogram is a simulation that includes LPM suppression, but not dielectric 
suppression, plus conventional transition radiation. The solid line includes LPM and di- 
electric suppression, along with conventional transition radiation. This was the 'standard' 
E-146 choice for simulation. 

Both suppression mechanisms are required to approach the data. However, there are 
still significant discrepancies between the LPM plus dielectric MC and the data. Below 800 
keV (350 keV) for 25 (8) GeV beams, the upturn in the data may be residual background, 
especially synchrotron radiation. The difference at higher photon energies is more complex, 
and will be discussed in the following subsection. 

Figure 17 shows the spectrum from 8 and 25 GeV electrons passing through 3% and 
6% X aluminum targets. The same three simulations are shown. Because aluminum has 
twice the Z of carbon, LPM suppression is considerably enhanced, with k^pu 8 MeV and 
800 keV at 25 and 8 GeV respectively. Because the density is similar to carbon, dielectric 
suppression is similar. Because dielectric suppression dominates for k < k cr , the curve for 
both suppressions is similar to that of carbon. The agreement between the data and the 
standard curve is better than with carbon. 

Figure 18 shows the spectrum from 8 and 25 GeV electrons passing through 3% and 6% 
X iron targets, with just the 'standard' MC. The general slope of the data matches the 
simulation, but the behavior at higher k for 25 GeV beams is quite different. 
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FIG. 16. Comparison of data from SLAC-E-146 with MC predictions for 200 keV to 500 MeV 
photons from 8 and 25 GeV electrons passing through 2% and 6% Xq carbon targets. The cross 
sections are given as dN / d(ln k) / Xq where N is the number of events per photon energy bin per 
incident electron, for (a) 2% Xq carbon and (b) 6% Xq carbon targets in 25 GeV electron beams, 
while (c) shows the 2% Xq carbon and (d) the 6% Xq carbon target in an 8 GeV beam. Three Monte 
Carlo curves are shown. The solid line includes LPM and dielectric suppression of bremsstrahlung, 
plus conventional transition radiation. Also shown are the Bethe-Heitler plus transition radiation 
MC (dashed line) and LPM suppression only plus transition radiation (dotted line). Adapted from 
Anthony et al. (1997). 
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FIG. 17. Comparison of data from SLAC-E-146 with Monte Carlo predictions for 200 keV to 500 
MeV photons from 8 and 25 GeV electrons passing through 3% and 6% X$ aluminum targets. The 
format and Monte Carlo curves are the same as in Fig. 16. Adapted from Anthony et al. (1997). 
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FIG. 18. SLAC-E-146 measurements and Monte Carlo predictions for 8 and 25 GeV electrons 
passing through a 3% and 6% Xq iron targets. The Monte Carlo curve is based on LPM and 
dielectric suppression, plus conventional transition radiation. Adapted from Anthony et al. (1997); 
Panel (c) is mislabelled as 6% Xq there. 
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FIG. 19. SLAC E-146 measurements and Monte Carlo for 3% Xq and 5% Xq uranium targets in 
8 and 25 GeV electron beams. The solid line shows the LPM and dielectric suppression, conven- 
tional transition radiation Monte Carlo prediction. The other lines include simulations based on 
calculations of transition radiation due to Pafomov (dashed line) and Ternovskii (dotted line) , Eq. 
(88), with x = 1- The flat solid line in panel (a) is a calculation based on Eq. (53). Adapted from 
Anthony et al. (1997). 



Figure 19 shows the bremsstrahlung spectra in uranium targets. Uranium is dense enough 
that LPM suppression is dominant, and the E-146 collaboration compared simulations with 
dielectric and LPM suppression, plus conventional transition radiation (TR), or the cal- 
culations of Ternovskii (1960) or Pafomov (1964) of transition radiation due to multiple 
scattering. 

The Pafomov curve jumps at k — k cr , around 800 keV (400 keV) at 25 (8) GeV, cor- 
responding from the switch from Eq. (89) to a numerical approximation for Eqs. (90) and 
(91). For k < k cr , Pafomov is considerably above the data and the conventional TR curve. 
Above the jump, the curve shows a reasonable trend, but the TR appears to be several times 
too high. 

The Ternovskii curve also jumps, at sk^/k 2 = 1, about 500 keV (below 200 keV) for the 25 
(8) GeV data. Below the jump, Ternovskii matches conventional TR. Above it, Ternovskii 
is quite far above the data. Moreover, it extends to too high an energy, above &lpm- 
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FIG. 20. SLAC E-146 data on 8 and 25 GeV electrons hitting a 0.7% X gold target. Shown are 
calculations by Blankenbecler and Drell (dashed line), and Shul'ga and Fomin (dot-dashed line). 
For comparison, the Migdal MC is shown as the usual solid line. Adapted from Anthony et al. 
(1997). 

The Ternovskii radiation could be reduced by lowering x below 1 in Eq. (88). However, a 
considerable adjustment would be required. 

Figure 20 shows the spectrum from 8 and 25 GeV electrons passing through a 0.7% 
X gold target. The solid histogram is the standard simulation. This target is especially 
interesting because, at 25 GeV, for k < 7 MeV, T < lf . Including LPM suppression, T < If 
for k < 3 MeV, and the target should interact as a single unit. The BH 1/k spectrum should 
be recovered, albeit at a reduced intensity. The dot-dashed line shows a calculation based 
on Eq. (53), in good agreement with the flattening observed in the data. For the E-146 
2% X lead (not shown here) and 3% X uranium 25 GeV data, Eq. (53) also applies, but 
only for a very limited range of k. However, as Fig. 19 shows, for 3% X uranium, Eq. (53) 
predicts dN/d(\nk) = 0.038 for 1.0 MeV < k < 3.7 MeV, considerably above the data. The 
lead data shows a similar discrepancy. Of course, these discrepancies may be because the 
equation is does not apply so close to T = If. 

Fig. 20 also shows the results of a calculation by Blankenbecler and Drell. Neither 
prediction has been corrected for multiple interactions in the target or detector resolution. 
Because neither calculation includes dielectric suppression or transition radiation, they fail 
for k < k cr . 

Fig. 21 shows the spectrum from 8 and 25 GeV electrons hitting a 0.07% Xo target. 
This target is only about 1.2 x (Xo/1720) thick, so that multiple scattering should produce 
little suppression, and the Bethe-Heitler simulation (dashed histogram) should be a good 
match to the data. Dielectric suppression should also be small, because the total phase 
shift in the target, h 2 ujpT/k 2 lf, is small for k > 500 keV. Transition radiation should also 
be strongly suppressed, by (T/lf) 2 . Without this suppression, transition radiation would 
completely dominate the data. As expected, the 25 GeV electron data matches the Bethe- 
Heitler spectrum. However, the 8 GeV data drops off for k < 2 MeV. For comparison, a 
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FIG. 21. Measurements and Monte Carlo for a 0.1% X gold target at (a) 25 GeV and (b) 8 GeV. 
The dashed line is the BH prediction (no suppression), while the solid line is a Monte Carlo which 
includes dielectric suppression, but not LPM suppression. Because the very thin target should 
exhibit little transition radiation, no transition radiation is included in these Monte Carlos. From 
Anthony et al. (1997). 

dielectric suppression only Monte Carlo simulation is shown. The data is considerably above 
this curve. 

Blankenbecler and Drell (Blankenbecler, 1997c) predict that, at 25 GeV, emission is 
suppressed by 10% compared to BH for k= 500 MeV, rising to 13% at k=100 MeV. At 
8 GeV, the difference is only a few percent. Unfortunately, these predictions differ from 
Bethe-Heitler by less than the experimental errors. 

For a target this thin, the collaboration has noted that the signal is very small and the 
potential backgrounds are large. Furthermore, the target thickness and overall normalization 
between the signal and simulations cannot be well determined. 

In all of these these plots, the Monte Carlo curves were normalized to the data by 
multiplication by a constant adjustment, chosen so the MC best matches the data above 20 
MeV (2 MeV) at 25 (8) GeV. The thresholds were chosen to avoid thin target corrections 
for T < Ifo, backgrounds and transition radiation; for the 0.7% Xo target, higher limits 
were chosen, 30 (10) MeV at 25 (8) GeV. Overall, the standard Monte Carlo curves had to 
be scaled up by an average of 5% (2a) to match the data. This discrepancy would likely 
disappear with an input cross section where the onset of suppression was more gradual 
around k ~ k LPM . 

The errors shown on the plots are statistical only. The E-146 collaboration has carefully 
studied the systematic errors on these measurements. The point-to-point systematic errors 
vary slowly with k and correspond to a 4.6% uncertainty for k > 5 MeV. Below 5 MeV, 
the systematic errors rise to 9%, because of increased uncertainties in photon energy cluster 
finding as Compton scattering takes over from showering as the dominant energy loss. The 
±3.5% systematic error on the normalization was determined separately. 
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FIG. 22. Bin by bin subtraction of the 3% Xq uranium data from the 5% Xq data by SLAC E-146, 
for 25 and 8 GeV beams. The solid line is the result when the same procedure was applied to their 
standard Monte Carlo. From Anthony et al. (1997). 



D. Bulk versus Surface Effects and Conclusions 

Although the data clearly demonstrates LPM and dielectric suppression, the thinnest 
targets show excess radiation over the bulk LPM predictions. This can easily be ascribed 
to surface radiation. As long as interference between the two target edges is negligible, the 
surface and bulk effects can be separated in a model independent way by subtracting spectra 
from targets made of the same material, but different thicknesses. 

Figure 22 shows the results of this subtraction for uranium. Because the subtraction 
exacerbates the effects of multiple interactions in the target, it is necessary to compare the 
subtracted spectra with similarly treated simulations. For all of the E-146 materials, the 
subtracted spectra and simulations agree more closely than the unsubtracted spectra, and, 
for most targets, the agreement is within the statistical and systematic errors. The notable 
exceptions are the 25 GeV iron and carbon data. 

The improved agreement indicates that the discrepancies in the unsubtracted spectra 
are due to surface effects. For the denser targets, this fits the expectations for transition 
radiation. However, for the lighter carbon and iron targets, the unsubtracted data indicates 
more suppression than expected. This implies that there is negative emission from the 
surfaces, an extremely unlikely result. 

The collaboration (Anthony et al, 1997) considered a few explanations for the discrep- 
ancies. For carbon, the increased suppression could be due to the crystalline structure of 
the pyrolitic graphite target. If the target density profile varied, then the higher density 
regions would show more suppression, and the aggregate would show somewhat higher LPM 
and dielectric suppression. For iron, magnetization of individual domains could produce 
some magnetic suppression. Unfortunately, the details depend on the (unknown) domain 
structure. 

The density profile could also vary if the target surfaces were oxidized. For the 2%X 
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tungsten target, the collaboration saw some evidence for this; thickness measurements based 
on weighing (giving results in g/cm 2 ) were lower than those obtained with calipers; an oxide 
layer could explain this. The anomalous normalization on the 2% tungsten target could also 
be explained by an oxide layer. For the lighter targets, a surface layer could explain the 
shapes of the curves, but the required thickness seems very unrealistic, and, in any case, the 
presence of such a layer is not supported by the normalization and thickness measurements. 

It is also possible that Migdal's theory may be inadequate for lighter targets. However, 
Zakharov (1998b) also found a very poor fit for the E-146 carbon data, despite a careful 
treatment of inelastic electron-electron interactions. In comparison, his fits to the other 
data were quite good. It is possible, although probably unlikely, that a different treatment, 
would find a better fit. It is also possible that a calculation taking into account the chemical 
structure of the pyrolitic graphite and using screening based on the actual electron density 
might better fit the data. 

The subtraction procedure could be modified to measure the transition radiation spec- 
trum from a single surface. However, the errors are slightly too large for the result to be 
interesting. 

Overall, the E-146 heavy target data is in good agreement with most of the recent 
treatments of LPM suppression. However, the shape of the spectra from the lightest targets 
remain a bit of a mystery. 

X. LPM EFFECT IN PLASMAS 

So far, we have considered particles moving through a medium. However, in a plasma, 
there is no separate incident particle, and each particle must be treated on an even footing. If 
the medium is sufficiently dense, there can be suppression even for non-relativistic particles. 
Then, the formation time tf is easier to use than the formation length. If tf is longer than 
the mean time between collisions t c = 1/T, T being the collision rate, then emission can be 
suppressed. 

One hadronic example of a plasma is a supernova; with the high temperature and density, 
t c is very short. For reactions with tf > t c , suppression may be present. One interesting 
reaction is the production of right handed neutrinos or axions through NN — > NNvT>, 
nn — > npeV e and NN — > NNa. If these hypothetical particles exist and are produced, they 
will carry energy away from the explosion, increasing the cooling rate. The measured cooling 
rate has been used to put limits on these particles (Raffelt and Seckel, 1991). 

For axions or neutrinos, tf = hfEE, where HE is the sum of the neutrino energies or the 
axion energy. When t c <C t/, then 'free' collisions are rare. Instead, the interacting nucleons 
are excited, This can be modelled by giving the nucleon an effective mass, as may be done 
with photons in dielectric suppression. Two complications come from the possibility of back 
reactions, and from the degeneracy of the incoming particles. The suppression is the ratio 
of the emission Q(T) to Q(T = 0): 



Here, n accounts for the multiplicity of the emitted particles, with n = 2 for axions and 
n = 4 for neutrino pairs; for axions this equation matches Eq. (42). Because T depends 




(104) 
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on the local density; finding the overall suppression factor for an entire supernova requires 
detailed modelling. However, Raffelt and Seckel estimate that S could be as small as 0.1. 
So, far fewer axions and neutrinos are emitted than if suppression were absent, and earlier 
limits that neglect this reduction are invalid (Raffelt and Seckel, 1991). 

The general solution for an non-equilibrium electromagnetic plasma has been elegantly 
formulated using non-equilibrium quantum field theory (Knoll and Voskresensky, 1995, 
1996). By careful classification of diagrams, and appropriate resumnation, Knoll and Voskre- 
sensky avoided infrared divergences, and reproduced both the classical (non-suppressed) and 
quasi-particle (low density) limits by appropriate choice of subsets of graphs. In a dense 
plasma, the quasi-free scattering approximation breaks down, and the reaction rate is re- 
duced by 

S -WTT- < 105 > 

This equation also matches Eq. (42), with jhup of the medium replaced by T, the relaxation 
rate of the source. 

A plasma like this might surround a quantum black hole. Before a quantum black hole 
explodes, it emits a huge flux of charged particles. This radiation forms a nearly thermal 
photosphere, consisting of electrons, positrons and photons (Heckler, 1995). The charged 
particle emission rate and consequent density is high enough that bremsstrahlung and pair 
production should be suppressed in this plasma. Since most of the emitted particles are 
hadrons, the black hole may also generate a dense color plasma, akin to a quark gluon 
plasma, with interactions suppressed as discussed in the following section. 

Similar plasma effects probably occurred during the big bang. Unfortunately, we do not 
know of any calculations involving it. 



XI. ELECTROMAGNETIC SHOWERS 

Although the LPM effect is best studied using single interactions in thin targets, most 
real world situations involve electromagnetic showers in thick targets. This includes nat- 
ural processes like cosmic ray air showers and neutrino induced electromagnetic showers, 
and showers in man-made detectors. Although modelling showers with suppression effects 
requires complex analytic calculations or Monte Carlo simulations, this section will present 
some simple calculations that show when suppression effects can be important. 

Suppression affects showers in several ways. Besides the obvious elongation when E > 
Elpm and the radiation length increases, the shower changes form, with the low energy 'fuzz' 
disappearing, and shower to shower fluctuations become much more important, because the 
number of interactions drops greatly. 



A. Natural Showers 

Many calculations have concerned high energy shower development in water or ice, an 
area relevant to high energy neutrino astronomy, or in air, an application motivated by 
extremely high energy (EHE) cosmic ray air showers. The latter case is complicated because 
the air density varies exponentially with altitude. 
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Water and ice showers are of interest for detecting very high energy astrophysical v e 
through resonant W production, v e e — > W — > ez/ e and similar reactions. This cross section 
rises rapidly at the W pole, corresponding to a v e initial energy of 6.4 PeV, an average of 2.1 
PeV of which goes to the outgoing electron. Since 2.1 PeV~ 7E LPM , LPM suppression is 
significant. As Fig. 2 shows, the electron radiation length is increased by about 35%. Since 
E ~ 4:E P , pair creation will also cause suppression for a range of k. 

At similar energies, the search for u T N — > tN — > ev T V e N events (Learned and Pakvasa, 
1995) in ice requires distinguishing the r track from the beginning of an electron shower; 
when suppression mechanisms reduce the number of low energy bremsstrahlung photons, 
this separation becomes more difficult. 

Because both applications involve very high energy showers, direct simulations have 
been limited by the available computer power; until recently at least partially analytic 
calculations were needed. These analytic methods owe much of their history to earlier 
analytic calculations of Bethe- Heftier shower development. 

One of the first analytic calculations (Pomanskii, 1970) showed that the penetrating 
power of electromagnetic showers rises as the LPM effect becomes important. At 10 19 -10 20 
eV in earth, electrons and photons become as penetrating as muons. Misaki (1990) used the 
matrix method to show that EM showers above 10 15 eV in water are elongated by the LPM 
effect. In the tails of the shower, at a given sampling depth, the LPM effect roughly triples 
the density. 

Unfortunately, while analytic calculations can determine the average shower shape, they 
have limited value for understanding shower to shower variations; for this, simulations are 
needed. To reduce the computing load, hybrid Monte Carlos are often used. A hybrid 
MC simulates the initial shower development; shower tails are added on based on a library 
of simulated complete showers. This reduces the computational requirements significantly. 
Stanev and collaborators (1982) measured shower elongation in water and lead due to the 
LPM effect using a hybrid MC. 

Konishi et al. (1991) used a hybrid MC to study fluctuations in a regime where sup- 
pression is very strong: 10 17 eV showers in lead. In contrast to BH showers, LPM showers 
showed large shower to shower variations. By eliminating soft bremsstrahlung, the LPM 
effect greatly reduces the number of interactions per radiation length, so the shower devel- 
opment depends on far fewer interactions, greatly increasing the shower to shower variation. 
The variation complicates energy measurement, especially for detectors with limited sam- 
pling. Misaki (1993), quantified the degree of fluctuation, measuring the distribution of the 
depth of shower maximum, the depth at which the number of shower particles is a maximum. 
For 10 17 eV showers in rock (Elpm — 77 TeV), he found that shower maximum occurs at 
18 ± 13 (FWHM) X for BH showers, compared with 152 ± 176 X for LPM showers. Not 
only is the maximum deeper, but its position varies much more. 

Not considered here is the increased angular spreading caused by suppression. Although 
it may be not affect the general shower shape, it is probably important in understanding 
radio waves produced by showers in ice (Zas, 1997), where the spectrum depends on the 
transverse separation between particles. 
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B. Cosmic Ray Air Showers 



Cosmic ray air showers occur when EHE cosmic rays hit the earths atmosphere and 
interacts to produce a cascade of particles (Sokolsky, Sommers and Dawson, 1992). Showers 
up to 3 x 10 20 eV have been seen (Bird et al, 1994). Two techniques are used to study the 
highest energy air showers. Large aperture telescopes, like the Flys Eye (Bird et al, 1994) 
observe the shower induced fluorescence of the N2 in air, measuring the shower development 
in the atmosphere. Ground based arrays of hundreds or thousands of small detectors, spaced 
up to a kilometer apart, observe the remnants of the shower that reach the ground (Takeda, 
1998) (Auger, 1996). 

Ground based detectors act as calorimeters with a single sampling layer, located behind 
a 28 Xq, or 15 hadronic interaction lengths, A thick atmospheric absorber. For a vertical 10 20 
eV shower, sea level is near shower maximum, so measurements are not too sensitive to the 
position of the first interaction. For obliquely incident showers, the detector is considerably 
behind shower maximum, and as a result, ground observations are subject to significant 
fluctuations depending on the position of the first interaction. Because of their limited 
ground coverage, most of their samples are many Moliere radii from the shower core, in the 
tails of the angular distribution. This heightens the sensitivity of the detector to the details 
of the initial interactions. 

On the other hand, air fluorescence detectors measure the whole shower profile, with the 
atmosphere acting as a fully active sampling medium. They are thus less sensitive to the 
details of the shower development. 

There has been considerable disagreement about whether the LPM effect is important 
in ultra-high energy showers. Capdevielle and Attallah (1992) found that LPM suppression 
had a significant effect on 10 19 to 10 20 eV proton induced air showers. However, Kalmykov, 
Ostapchenko and Pavlov (1995) found a much smaller change; for a 10 20 eV shower, they 
found a 5% decrease in the number of electrons at shower maximum, and a 15 ± 2 g/cm 2 
downward shift in the position of shower maximum. Although this shift is smaller than the 
measurement resolution, it is systematic, and important in studies of cosmic ray composition. 
Kasahara (1996) found that the LPM effect reduces the number of particles reaching the 
ground by about 10% for a 5 x 10 20 eV proton shower. None of these authors considered 
the effect of other suppression mechanisms. 

Although a complete Monte Carlo simulation is required to quantify the effect of sup- 
pression, simple calculations can illustrate some qualitative features of shower development 
(Klein, 1997). The calculations depend strongly on the identity of the incoming particles, 
with protons (or neutrons) the most popular, although photons or heavier nuclei cannot be 
excluded. For heavy ions, suppression is greatly reduced because of the lower per-nucleon 
energy. Because this model is very simple, the possibility of photons pair converting in 
the earth's magnetic field will be neglected (Kasahara, 1996)(Stanev and Vankov, 1997). 
The probability depends on the photon energy and the angle between the photon and the 
earths magnetic field. For primary photons with k > 10 20 eV, the probability is large, partly 
because the magnetic field extends to much higher altitudes than the atmosphere. 

Because the atmospheric pressure, l/X and 1/Elpm decrease exponentially with height, 
it is convenient to use an exponential depth variable such as column density. Column density 
is measured in g/cm 2 , with ground level, A at 1030 g/cm 2 . Then, E LPM = (A /A) 234 
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PeV, where A is the depth in column density. Similarly, y die = 1.3 x 10~ 6 (A/A ) and 
E p = ^JAq/A 42 PeV. Because the atmosphere is much cooler at high altitudes, these 
numbers underestimate suppression by about 25%. They also neglect the fact that air is 
composed of diatomic N 2 and 2 molecules; when an electromagnetic interaction involves 
one atom of the molecule, the other atom will introduce additional multiple scattering over 
that expected from a monoatomic gas. This could be a significant effect, but it has yet to 
be studied. 

Incoming photons react by pair production, while protons interact hadronically. A central 
hadronic collision will produce a shower of several hundred pions; the neutral pions will decay 
to photons. The highest energy 7r° will have a rapidity near to the incoming proton, and 
their decay photons will have energies up to 2 x 10 19 eV. Many diffractive processes, such as 
A production can also produce photons with similar energies. Overall, photons from central 
interactions will have an average energy of about 2 x 10 17 eV. Above ~ 10 18 eV, the finite 
7T° lifetime becomes important (Kasahara, 1996); a 2 x 10 19 eV ir° travels about 5 km before 
decaying and is likely to interact hadronically instead of decaying. This will reduce the 
number of high energy photons, reducing the significance of LPM suppression in the earliest 
part of the shower development. 

Because both E^pm and the average particle energy, E, decrease with depth, suppression 
mechanisms can actually become stronger as a shower moves deeper in the atmosphere. 
Fig. 23 compares E LPM with E for an idealized Bethe-Heitler shower from a 3 x 10 20 eV 
photon. In each successive radiation length, there are twice as many particles with half the 
energy. The dashed line shows a similar cascade, from a 2 x 10 19 eV photon starting at 1A. 
The electromagnetic interactions at a given depth are determined by the ratio E/Elpm- For 
the photon shower, suppression is largest around 75 g/cm 2 , where E ~ AQE^pm- Electron 
dE/dx is reduced about 80%, and pair production cross section is reduced by 60%; X has 
more than doubled. For hadronic showers, the effect is smaller, and, of course, these high 
energy photons are only a small portion of the total shower. On the other hand, since 
hadronic interactions are only partially inelastic, the proton may carry a significant fraction 
of its momentum deeper in the atmosphere, where suppression is larger. In general, at least 
part of the shower will be suppressed. Because of the large shower to shower variations, it 
is difficult to give more quantitative estimates. 

This simple model underestimates the importance of suppression. When suppression 
slows shower development, E of the remainder of the shower will increase, increasing the 
suppression in the next radiation length. However, the model shows that suppression is very 
important in photon showers, and in at least parts of proton initiated showers. However, 
suppression is clearly less significant for proton showers than predicted by Capdevielle and 
Attallah (1992). These estimates appear consistent with Kalmykov, Ostapchenko and Pavlov 
(1995) and Kasahara (1996). 

Fluctuations are also very important in air showers. Because the cosmic ray energy 
spectrum falls as dN/dE 1/E 3 , it is important to understand the tails of the energy 
resolution distribution; without accurate simulations, showers whose energy is overestimated 
can skew the measured spectrum. Suppression mechanisms exacerbate the fluctuations, by 
greatly decreasing the number of interactions in the initial stage of the shower. For example, 
at sea level, a 10 10 eV electron emits an average of 14 bremsstrahlung photons per X , while 
a 10 17 eV electron emits only 3 photons. Pair creation is similarly affected; pairs become 
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FIG. 23. E LPM (solid line), E for a 3 x 10 20 eV photon shower (long dashes) and E for a 2 x 10 19 
eV photon created at 1A (short dashes). Es is based on an idealized Bethe-Heitler shower where the 
number of particles doubles each radiation length. At altitudes where E > Elpm, suppression is 
important; the ratio of the two energies determines the degree of suppression. Here, a temperature 
correction has been added to Elpm- 
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more and more asymmetric. The result is that shower to shower fluctuations are much 
larger. 

Suppression also greatly reduces the number of particles in the early stages of a shower, 
eliminating most of the low energy 'fuzz.' The LPM effect, dielectric suppression and pair 
conversion suppression all contribute to this reduction. Although the LPM effect covers 
the widest range of energies, the other effects are stronger in the limited range in which 
they operate. For E > E p , dielectric suppression reduces the number of bremsstrahlung 
photons with k < 331 MeV by two orders of magnitude. None of the current simulation 
efforts include these other mechanisms. Even relatively late in the shower, there will be a 
reduction in the number of low energy particles. For example, where E ~ 10 13 eV, the LPM 
effect suppresses photons below ~500 MeV. Although a low energy cutoff does not affect 
the overall shower development, it can significantly reduce the number of particles early in 
the shower. This reduction may affect the shower profile observed by future air fluorescence 
detectors. 

Because Elpm drops as altitude decreases, suppression greatly increases the chance of 
photons penetrating deep in the atmosphere, even if the average interaction depth isn't too 
different. For trajectories where the probability of pair conversion in the magnetic field is 
small, Bethe-Heitler and LPM predict similar average conversion depths for a 3 x 10 20 eV 
shower, 114 and 122 g/cm 2 respectively. However, Migdal predicts that the photon has a 
7% chance of surviving to a depth of 6X , while Bethe and Heitler predict only a 0.25% 
survival probability at the same depth. These occasional deep interactions produce showers 
quite different from photons that convert higher in the atmosphere, and much more energy 
will reach the ground. Although the magnitude is smaller, similar effects may be seen for 
the most energetic photons from proton interactions. 

There have also been searches for neutrino interactions with air deep in the atmosphere, 
producing nearly horizontal air showers. Because the neutrino interaction cross section rises 
with energy, this becomes especially attractive for very high energy neutrinos, E ~ 10 15 — 10 21 
eV (Capelle et al, 1998). Since these neutrino interactions occur at low altitudes, where 
Elpm is a few hundred PeV, suppression is extremely important, especially at the higher 
energies. For example, a 10 20 eV electron at sea level has a dE/dx a factor of 10 lower than 
the Bethe-Heitler prediction. 

Although suppression mechanisms may be of little importance to current cosmic ray stud- 
ies, the next generation of experiments will reach both higher energies and higher statistics 
(Auger, 1996, Krizmanic, Ormes and Streitmatter, 1998). Although careful, detailed sim- 
ulations are needed to reach definite conclusions, it appears that suppression mechanisms 
will affect the data collected. 

C. Showers in Detectors 

Changes in shower energy deposition may also be important in the next generation of 
high energy and astrophysics detectors. This review will consider two examples, with a few 
representative numbers. 

Direct cosmic ray composition measurements depend on measuring the cosmic ray charge 
and energy; the latter in a calorimeter. Because these experiments must be done above much 
of the atmosphere, either in balloon flights or space based experiments, calorimeter thickness 
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is limited by the allowable weight. For example, the JACEE collaboration has flown ~ 6X0 
thick lead + emulsion/X-ray film calorimeters on long duration balloon flights (Parnell et 
al, 1989). Protons with energies up to 500 TeV were observed. Since a proton can be 
diffractively excited to a A particle, which can decay pn° — > ^77, photons with energies 
up to about 7% of the proton energy, 35 TeV, can be produced. Figure 2 shows that the 
location of the photon conversion point is not affected since 35 TeV ~ 8Elpm- However, 
the produced electrons and positrons will exhibit reduced energy loss, and overall energy 
leakage out of the back of the calorimeter may increase. The next generation of space 
based experiments will detect protons with energies well beyond 1 PeV, with a consequently 
increased effect (Parnell et al, 1989). 

At the Large Hadron Collider, now under construction at CERN, electrons with energies 
up to about 1 TeV will be produced. Electron identification will be aided by a pre-shower 
radiator, which will separate electrons from hadrons by measuring the energy radiated in the 
first 2 or 3 X of a lead calorimeter (Aspell et al, 1996). Because suppression reduces the 
number of photons emitted by the electron, electrons will behave more like heavier particles 
such as muons and pions. Without suppression, the number of photons emitted by a 1 
TeV electron in lead is infrared divergent: iV 7 ~ T / Xq In (E/k m i n ). For a detection cutoff 
Emin = 1 MeV, this is 14 photons per radiation length. Suppression reduces this to an 
average of 3 photons per radiation length. The small average number of photons increases 
the probability that the early stages of an electron shower will be indistinguishable from a 
muon or pion interaction. 

XII. QCD ANALOGS 

QCD analogs of the LPM effect involve quarks and gluons moving through a strongly 
coupling nuclear medium. Because the nuclear medium is so dense, a quark or gluon moving 
through the medium typically undergoes many soft color exchange interactions, analogous 
to multiple scattering. Because a nucleus is small, the formation length is often larger than 
the nucleus, so the nucleus acts as a single radiator. 

The literature on this subject is voluminous and this is not the place for a detailed review. 
Instead, we will discuss a few examples, to illustrate the basic features of the interactions, 
and also to emphasize the both the similarities and the differences between the QED and 
QCD cases. Because QCD couples strongly and is non-Abelian, the theoretical calculations 
are much less robust, and the experimental data is much harder to interpret. 

A. Hadron Level Calculations 

In hadronic interactions, particle masses play a larger role than they do in electromagnetic 
interactions. Suppression effects in hadronic interactions can be seen easily by considering 
the photoproduction of vector mesons, for example — > <pN . This reaction is similar to 
7 — > e + e~, with a quark pair produced instead of a lepton pair. One difference is that the 
interaction between the quark pair and the nucleus is hadronic, rather than electromagnetic. 

The formation length for vector meson production is the same as for pair production, 
2hk / MyC 3 . For a 100 GeV photon producing a 0, lf = 40 fm, far larger than the nuclear 
diameter. So, the photon/vector meson interacts with the nucleus as if it is a single particle, 
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with the interactions determined by Eq. (53). Unfortunately, QCD does not specify f{6), 
but, in principle, this equation could be used to explore how the cross section depends on 
the nuclear thickness. 

Another example where suppression effects can be important is low mass Drell-Yan 
dilepton production, qq — > l + l~ . In the rest frame of the target nucleus, the incoming quark 
or antiquark can be thought of as pair producing the leptons. The momentum transfer from 
the target quark or antiquark can be small, so the formation zone can extend far outside 
the nucleus. 



B. Quark Level Calculations 

More fundamental calculations involve quarks and gluons moving through a possibly 
very hot nucleus. One focus of recent calculations is the search for the quark gluon plasma; 
a fast quark in a hadron gas (normal nucleus) may interact differently than one in a quark 
gluon plasma in which protons and neutrons are replaced by a sea of quarks and gluons. 
The sea of quarks and gluons, will have a longer screening length and a larger cross section, 
so quarks and gluons should lose energy considerably faster than in medium consisting of 
confined quarks and gluons. Measurements of high energy jets (or hadrons) energy loss has 
been proposed as a tool for detecting the quark gluon plasma (Wang, Huang and Sarcevic, 
1996). 

It is instructive to consider bremsstrahlung energy loss (dE/dx) by a quark with energy 
E emitting gluons of energy A; in a nuclear medium. Including the quark and gluon masses 
m q and m g (S0rensen, 1992) 

/ 2hEk{E - k) rintt 

f° myk 2 + myE(E-k)' V } 

Unfortunately, the quark and gluon masses critically affect Ijq. Besides bremsstrahlung, 
gluons can 'pair create' quark anti-quark pairs, a close analog of pair production. Because 
of the possibility of experimentally measuring jet or hadron energy loss discussed above, 
most calculations have focused on bremsstrahlung. 

S0rsensen (1992) repeated the Landau and Pomeranchuk semi-classical derivation, and 
also adapted Migdal's formulae to quarks and gluons. Although the results depend on m q 
and m g , and hence have large systematic errors, he did show that, for reasonable mass 
choices, suppression was important. This mass problem can be avoided because, for QCD, 
the gluon emission angle (equivalent to 9 7 is large, and Eq. (37) can be written as 

i f o = |r (107) 

where k± is the perpendicular energy of the emitted gluon. This avoids any explicit mass 
dependence. 

Brodsky and Hoyer (1993) used quantum mechanical arguments to find the energy loss 
of a parton travelling through nuclear matter. They decomposed the incident hadron into 
a spectrum of Fock states of varying masses, thereby avoiding the parton mass uncertainty. 
The resulting energy loss is limited by AE / E < kA 1 ^ 3 ' jx\s where k ~ 0.5 GeV 2 , x\ is the 
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fractional energy of the produced parton and yfs is the center of mass energy. Here, Xiy/s is 
analogous to E. The target thickness is given by the nuclear radius, A 1 / 3 . However, surface 
effects are not included in the calculation. Numerically, this energy loss should be negligible 
for reasonably high energy hadrons. 

Wang, Gyulassy and Plumer (1995) did a detailed calculation of energy loss in a quark 
gluon plasma. They found that the radiation, while obeying the bound listed above, was 
very sensitive to the color screening distance in the plasma. The strength of the suppression 
depends on the ratio of the mean free path A to the formation time tf, as with the QED 
plasma. When tf < A/c, 



where a s is the strong force coupling constant, (qj_) is the square of the average momentum 
transfer from a single scattering, which should be proportional to fi 2 , with /i the color 
screening mass. For quarks, the color factors C2 = 4/3 and r 2 = 9/8, while for gluons 
C2 = 3 and r 2 = 9/8. So, the energy loss only grows logarithmically with particle energy; a 
typical value is dE/dx ~ 3.6 GeV/fm. In the opposite limit, tf > A/c, dE/dx ~ E, as with 
Bethe Heitler. 

More recent calculations have emphasized the finite size of nuclei. Many of these works 
have also considered the QED case, useful as a check of their results. Zakharov (1997) applies 
the techniques of Sec. V to finite thicknesses of both nuclear matter and the quark gluon 
plasma. For a quark travelling through a nuclear medium, the energy loss is proportional 
to the distance travelled, AE/E ~ T/10 fm, considerably higher than Brodsky and Hoyer. 
Zakharov found that surface effects are large, so LPM suppression plays a limited role. This 
probably explains the disagreement with Brodsky and Hoyer, who neglected surface effects. 

For high energy quarks (Z/o > Ra) created inside the nucleus, AE ~ T 2 ; this dependence 
comes about because the quark takes some time to develop its gluon field; this is analogous to 
the lower dE/dx observed for newly created e + e~ pairs. For quarks created at lower energies, 
where If < Ra, the energy loss is proportional to the distance travelled, AE/E ~ T/10 fm. 

R. Baier and collaborators (1995) considered a fast quark or gluon propagating through 
nuclear matter. In an infinite medium, the soft gluon spectrum matches their soft photon 
QED result (R. Baier et al. 1996), with dE/dx ~ \[E. This differs from the logarithmic 
energy dependence found by Wang, Gyulassy and Plumer. The difference is that this work 
included additional diagrams, such as the non-Abelian 3-gluon vertices; they note that these 
are the dominant source of emission. A later, expanded collaboration (R. Baier et al., 1997) 
considered a finite size quark gluon plasma and found the same energy dependencies as 
earlier, along with the same R\ dependence found by Zakharov. However, for lower energy 
emission, they find AE/E ~ T/4 fm, much higher than Brodsky and Hoyer or Zakharov. 
This is because Zakharov includes additional diagrams in his calculation, and also treats 
hard gluons differently (R. Baier et al., 1998b). 

The collaboration has recently studied energy loss in a longitudinally expanding quark 
gluon plasma (R. Baier et al, 1998a). This is the first attempt to model a realistic, time 
varying temperature and density distribution that could occur when two heavy nuclei collide. 

While the different calculations agree in many ways, there is still some significant dis- 
agreement. They generally agree about the appropriate energy scaling in the Bethe-Heitler 




(108) 



67 



(no suppression) and strong suppression regimes, and also show a good correspondence with 
the QED calculations. The disagreement is over where these two regimes apply. Some of this 
stems from differing treatment of surface terms. Some of it may stem from the details of the 
initial state used. Whatever the cause, the numerical results vary greatly. Unfortunately, 
because of the difficulty of clear experimental tests, it may be some time before data can 
choose the best result. 

XIII. SUPPRESSION IN E+E COLLISIONS 

Future high energy electron-positron colliders will collide extremely dense beams of elec- 
trons and positrons. Besides the desired hadronic interactions, large numbers of photons 
will be created. Photon emission before the hadronic interaction lowers the average collision 
energy. Also, some of the photons will pair convert, introducing a charged particle back- 
ground into the detector. Photons can be produced by coherent beamstrahlung, where a 
single particle interacts with the other beam as a whole, or by incoherent bremsstrahlung. 
Because of the density of the plasma, multiple scattering, Compton scattering and magnetic 
suppression can suppress some of these interactions. 

Although suppression occurs for similar reasons as in bulk matter, many of the details are 
different because of the very different environment. Both bremsstrahlung and pair creation 
are modified because bare charges are involved; q min is governed by the size of the beams. 
For both, the magnetic fields produced by the opposing beam can have strong effects, and 
both bremsstrahlung and pair creation are significantly suppressed. (Baier and Katkov, 
1973)(Katkov and Strakhovenko, 1977), as is beamstrahlung. 

In contrast, because of the low density, LPM suppression is small in currently envisioned 
machines. Dielectric suppression transfers over rather directly, after adjusting for the dif- 
fering electron densities. It is in principle strong enough to cause significant suppression. 
However, because the beams are much shorter than lf , dielectric suppression is reduced, 
and of only marginal significance (Chen and Klein, 1993). 

XIV. OPEN PROBLEMS AND FUTURE POSSIBILITIES 

Theoretically, there has been much progress in the past few years, and several new 
approaches to LPM suppression have appeared. The problem of radiation from slabs, sets 
of foils and the like have all been considered for a wide range of kinematic factors. Despite 
this, these calculations have some limitations. Most only consider multiple scattering, and 
none cover the full range of suppression mechanisms. None cover pair production; this may 
no longer be a trivial extension. 

Finally, the high energy regime when bremsstrahlung and pair creation suppress each 
other needs to be explored. Even if QED tests are not practical, this regime may be im- 
portant for understanding the QCD case. More generally, none of the calculations consider 
higher order diagrams. 

Even in the absence of a unified approach, a better treatment of dielectric suppression 
is of interest. To treat multi-GeV bremsstrahlung with classical electromagnetism concepts 
like the dielectric constant of a medium is rather anachronistic. The dielectric constant 
comes from photons forward Compton scattering off the electrons in the medium; it would 
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be very nice to see a calculation of dielectric suppression that uses Compton scattering as a 
starting point; such a calculation might yield some interesting surprises. 

All of these recent results need to be made generally accessible. A public program library 
would let experimenters compare the different approaches under a variety of conditions. 

Many implications of the LPM effect are poorly appreciated, and there is a need to 
learn more about them. One 'lesson' from SLAC-E-146 is that there is very little data 
on small y bremsstrahlung, and consequently little theoretical attention to this area. Of 
particular concern are radiative corrections to electromagnetic effects. Besides the effect 
that suppressed bremsstrahlung will have on elastic scattering, little is know about when 
higher order terms become important. 

Experimentally, there is a need to accurately measure suppression at higher energies, and 
to look at new phenomena. Cosmic ray based experiments have poor statistics, and current 
accelerators only reach the semi-classical (k <C E) regime. Ideally, a study of quantum 
bremsstrahlung would use electrons with E > E LPM ks 2.5 TeV for gold. An experiment 
at Fermilab (Jones et al, 1993) with 250 GeV electrons would be a significant step forward, 
but would not reach this goal. Unfortunately, a TeV electron beamline seems a long way 
off. However, because a relatively low flux is sufficient, it might be possible to use electrons 
and or photons produced in 14 TeV pp collisions at LHC to study bremsstrahlung or pair 
production. Suppression effects are large enough that they may also be observable in the 
LHC general purpose detectors. This is particularly likely in apparatus such as pre-shower 
radiators that look at early shower development. 

Any future test of LPM or dielectric suppression must at least match E-146 in statis- 
tics, backgrounds and control of systematics. The statistics should be simple, but the 
backgrounds and systematics will require some effort. It will be necessary to minimize the 
magnetic fields in the bending magnets, photon angular acceptance, and backgrounds due 
to the beam transport system. With a higher energy beam, it should also be possible to 
reduce the systematic errors, by avoiding the calorimetric 'transition region' between en- 
ergy deposition by Compton scattering and by showering. It would also be very useful to 
collect data with a wider range of target thicknesses, and also more low Z targets. With 
an experiment accurate to 1-2%, it should be possible to compare the data with the recent 
theoretical predictions, assuming that the problems of multiple emission by a single electron 
can be handled properly, either theoretically or experimentally. While experiments with 
higher accuracy than E-146 are very desirable, improving the systematic errors by a large 
factor over E-146 would require extreme attention to detail. 

Any new experiment should also investigate other phenomena, such as magnetic suppres- 
sion. With E > 200GeV, ys > Udie and magnetic suppression should be directly observable. 
Studies in a constant magnetic field are probably the most interesting, but the effect of 
randomly ordered domains might also be observable. Synchrotron radiation could be sep- 
arated from bremsstrahlung by measuring targets of different density in a fixed magnetic 
field. A similar apparatus could measure emission from targets consisting of stacks of foils, 
to investigate transition radiation due to multiple scattering. It is desirable to try to push 
bremsstrahlung measurements to the smallest y possible. As Eq. (102) demonstrates, this 
is a new regime, and some surprises may await us. One tangible goal would be to reach the 
point where higher order corrections cause current theories to fail. 

There is currently considerable theoretical effort studying LPM-like phenomena in QCD. 
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With the coming data on high energy heavy ion collisions from RHIC, this work will probably 
continue, with an increasing emphasis on ways to differentiate between suppression in a 
hadron gas and that of a quark gluon plasma. Also, the current calculations are limited by 
their neglect of higher order terms, which are likely to be very important. Experimentally, 
there is a clear need is for a convincing demonstration of suppression involving QCD. 

One place where more theoretical work is needed is to investigate suppression in astro- 
physical phenomena. One down-to-earth need is a detailed study of suppression in EHE 
cosmic ray air showers. Other topics, less tied to specific experimental techniques, include 
the study of suppression in astrophysical regions of extreme density, temperature, and/or 
magnetic field. This review has discussed a few such areas, but there are many more. 

Between the new theoretical approaches and the expanding range of applications in QCD, 
plasmas and astrophysics, suppression mechanisms seem destined to be a growing area of 
study over the coming years, 
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